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Abstract. Rapid progress has been made during recent years in the 
understanding of the flavour oscillations that occur as neutrinos traverse through 
supernova. The previous paradigm has given way and it is now clear that the 
neutrino signals we shall receive from future Galactic supernovae will allow us 
both to peer inside these extraordinary cosmic events and to probe some of the 
fundamental properties of these elusive particles. In this review we aim to distill 
the progress that has been made focusing upon the effects of the dynamic density 
profile and the emergence of collective flavour oscillations due to neutrino self- 
interactions. 
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1. Introduction 

The detection of neutrinos from the Sun (Davis et al, 1984, 1968) and Supernova 
1987A (Alekseev et al, 1987; Bionta et al, 1987; Hirata et al, 1987) have laid the 
foundation for a new field, neutrino astronomy. These discoveries not only opened a 
new window upon the Universe through the observation of cosmic neutrinos but also 
established a new way to probe fundamental properties of these elusive particles in the 
grand "laboratories" of stars or even the entire Universe. One such experiment was 
the Sudbury Neutrino Observatory (SNO) (Ahmad et al, 2002) which showed that 
neutrinos other than ly^, the only neutrino species produced in the Sun, were present 
in the solar neutrino flux at Earth. This transformation to other flavours and the 
concomitant deficit in Ve — the solar neutrino problem (Bahcall and Sears, 1972) — 
can be elegantly explained by the Mikheyev-Smirnov-Wolfenstein (MSW) mechanism 
(Mikheyev and Smirnov, 1985; Wolfenstein, 1978). If neutrinos have masses and the 
flavour states are not the mass eigenstates, then the electron neutrinos emitted in the 
high density at the core of the Sun experience resonant flavour conversion as they 
propagate to the vacuum thereby altering the composition of the flux. SNO, other 
solar neutrino experiments, atmospheric neutrino experiments and also terrestrial 
experiments have now compiled compelling evidence that indeed neutrinos of different 
flavours can mutate into each other and this phenomenon is known as "neutrino flavour 
transformation" or "neutrino (flavour) oscillation" (see, e.g. Amsler et al, 2008 for a 
review) . 

But the focus of this review is core-collapse supernova, which we will refer to 
simply as supernova from now on. The explosion marks the death of a star with 
mass greater than ~ 8Mq (see, e.g. Woosley et al, 2002 for a review) initiated as the 
iron core of the star collapses when it is incapable of generating sufficient supporting 
thermal energy through further nuclear fusion (Weaver et al, 1978). A proto neutron 
star (PNS) (Baade and Zwicky, 1934) forms at the centre of the supernova if the 
collapse can be halted after supra nuclear density is reached. Most of the gravitational 
binding energy (~ 3 x 10^'^ erg) of the PNS is released in the form of neutrinos of all 
species with energies ~ lOMeV over a period of 10 seconds or so. Supernovae are 
such luminous sources of neutrinos that a supernova located 10 kpc from Earth would 
flood a neutrino detector here with a flux of ~ 10^^ cm^^s^-'^. Over the duration of 
the burst this flux is more or less the same as that from the Sun but with much larger 
average energy. With current neutrino detector technology we should observe at least 
several thousands from the next supernova in our Galaxy (Ahrcns et al, 2002; Beacom 
et al, 2002; Cadonati et al, 2002; Dornic et al, 2008; Ikeda et al, 2007; Sharp et al, 
2002). 

Because supernova energetics are dominated by neutrinos, it was immediately 
realized after the discovery of the MSW mechanism that any neutrino flavour 
transformation in supernovae could have important effects (Fuller et al, 1987). There 
are many flavour transformation scenarios that could occur in the supernovae but 
for this review we shall consider just the flavour transformation of active neutrinos 
outside the neutrino sphere where neutrinos start free streaming. In this regard the 
flavour transformation is similar to that in the Sun but there are three key aspects 
of neutrino flavour transformation in a supernova that also make it very different 
from the solar environment: (a) the matter density in a supernova is much higher 
{p ^ lO^^gcm^'^ at the neutrino sphere) so that oscillations can be affected by both 
neutrino mass- squared differences; (b) the density profile of a supernova is rich with 
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prominent features which evolve over the duration of the neutrino emission; and (c) 
the neutrino flux near the neutrino sphere can be so large that neutrino self-interaction 
must be taken into account. The vast majority of the neutrinos emitted by the PNS 
propagate through the overlying material unimpeded; only a small percentage will be 
absorbed by nucleons/nuclei. But due to neutrino oscillations what is detected here 
on Earth is not the same as that emitted by the PNS. The neutrino evolves as it 
propagates and for practical purposes the evolution to Earth is solved in a sequence 
of steps: first, we calculate the state at the surface of the supernova; second, we 
add on the effect of propagating through the vacuum to Earth; and third, an "Earth 
matter" effect (Dighe and Smirnov, 2000) is sometimes included to account for the 
possibility that the neutrinos may have to pass through the Earth before reaching 
a detector. Of the three steps in calculating the neutrino signal at Earth it is the 
first step, the calculation to the surface of the supernova, that is by far the most 
difficult and which forms the focus of this review. We should mention here that at the 
present time neutrino oscillations are not part of supernova simulations or are treated 
only approximately. So to date neutrinos are usually sent through the supernova in a 
post-processing stage, i.e. the supernova is an input. 

1.1. MS W flavour transformation with dynamic density profiles 

The first studies of the electron dominated/pure MSW effect in supernovae employed 
purely static density profiles. Temporal evolution of the profile was ignored and any 
time dependence of the neutrino signal was limited to the variation of the neutrino 
spectra emitted by the core. But in recent years this has changed and the density 
profiles used now have become increasingly sophisticated. 

The first, generic, evolutionary feature of any core-collapse supernova is the 
forward shock. That a shock is formed in the supernova was discovered in the first 
simulations by Colgate and Johnson (1960) and the feature is generated by the rebound 
of the PNS during the implosion (Colgate et al, 1961). Initially it was thought that the 
shock would then proceed to shed the mantle of the supernova, the so-called prompt 
explosion mechanism, but as the physics input into the simulations became more 
accurate it was eventually realised that actually the prompt explosion mechanism fails 
(Hillebrandt and Mueller, 1981; van Riper, 1982). The outward motion of the forward 
shock stalls at a distance ~ 200 km from the PNS chiefiy due to nuclear dissociation 
(Arnett, 1982; Mazurek, 1982). If the star is to explode the shock must be revived from 
its stalled position and this can only occur by heating the material behind the shock. 
Due to the prodigious amounts of energy they convey, the neutrinos streaming out 
from the PNS were invoked by Bethe and Wilson (1985) as the transporters of energy 
from the hot core to the region behind the shock. This current paradigm of a stalled 
shock that is later revived is termed the delayed explosion mechanism. Revival of a 
shock due to neutrino heating does appear to work in the limited case of stars in the 
mass range 8-11 Mq which have 0-Ne-Mg cores (Dessart et al, 2006; Kitaura et al, 
2006) . The explosion energies of the simulations are smaller than the typical observed 
energies of supernova by around an order of magnitude but such weak supernovae 
do seem to match the energetics of some faint, ^^Ni deficient. Type II P supernovae 
(Chugai and Utrobin, 2000; Pastorello et al, 2007). But for more massive stars where 
the core is composed of iron all spherically symmetric simulations of the supernova 
using Boltzmann neutrino transport do not explode (Hix et al, 2003; Liebendorfer 
et al, 2001; Mezzacappa et al, 2001; Rampp and Janka, 2000; Thompson et al, 2003). 
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Any simulation that does explode uses "gray" neutrino transport (Mezzacappa, 2005) . 
While less satisfying than an explosion based upon first principles such simulations are, 
nonetheless, useful because we can learn from them which features we should expect to 
find when truly successful supernova simulations are performed. Schirato and Fuller 
(2002) — and later Lunardini and Smirnov (2003), Takahashi et al (2003), and Fogli 
et al (2003) — were the first to consider the effect of the dynamic profile upon the 
neutrinos. They demonstrated that the forward shock would race out through the 
mantle of the star and reach the H resonance region — a position we shall clarify 
in the next section — within the first couple of seconds. Upon arrival its presence 
in this region alters the adiabaticity of the neutrinos as they propagate through the 
supernova and thus the flavour content of the neutrino signal at Earth. This discovery 
that the evolution of the density profile alters the flavour composition of the neutrino 
signal is the reason that one now needs to use time dependent density profiles in order 
to make accurate predictions. 

After Schirato and Fuller (2002) additional features in the profile were also found 
to infiuence the neutrinos. The energy deposition does not switch off once the forward 
shock is revived and has begun moving outwards again. The material close to the PNS 
continues to be heated and the effect is to create a wind with a velocity that increases 
with distance from the core (Duncan et al, 1986). This wind pushes against the slower 
moving material in front of it creating a bubble or cavity. If the post shock-revival 
energy deposition is even more vigorous the wind strength grows to the point where 
its velocity may supersede the local sound speed. This leads to the formation of a 
new shock in the profile that faces the PNS rather than the exterior (Burrows et al, 
1995; Janka and Mueller, 1995), i.e. it is "reversed". Compared to the forward shock, 
the reverse shock is much more skittish. The reverse shock can be so strong that it 
penetrates to densities below those in front of the forward shock and then later it can 
diminish in size as the wind abates and even stall and return to the core (Arcones 
et al, 2007; Kneller et al, 2008; Tomas et al, 2004). The effect upon the neutrinos of 
the presence of the reverse shock in the the profile was first considered by Tomas et al 
(2004) . Unfortunately these authors unintentionally omitted the consequences of the 
correlation of the wavefunction between the multiple H resonances now present in the 
profile which can lead to phase effects. That multiple resonances can lead to phase 
effects had been discussed many years prior by Kuo and Pantaleone (1989b) but it 
was only recently that Fogli et al (2003), re-examining the profiles of Schirato and 
Fuller (2002), discovered the phase effects in supernova neutrinos. Further examples 
showing their presence were presented by Kneller and McLaughlin (2006) and then 
phase effects were discussed extensively by Dasgupta and Dighe (2007). 

But the failure of the best one-dimensional (ID), first-principles simulations 
indicates that supernova must be a multi-dimensional phenomena requiring rotation, 
convection, magnetic fields and/or other multi-dimensional physics. If that is the 
case then one would expect any explosion to be aspherical with a density profile that 
varies with the line of sight. Observational evidence that supernova are aspherical has 
abounded for years: asphericity has been observed or is implied from observations of 
Supernova 1987A (Dotani et al, 1987; Erickson et al, 1988; Li et al, 1993; Matz 
et al, 1988; Sunyaev et al, 1987), from the Vela SNR (Aschenbach et al, 1995), 
from Cassiopeia A (Hughes et al, 2000; Hwang et al, 2004) and from many other 
supernova remnants (e.g. Katsuda et al, 2008; Leonard et al, 2006; Wang et al, 
2001). And of course pulsars have been observed with astounding velocities with 
some approaching ~ 1500 kms~^ (Chatterjee et al, 2005; Cordes et al, 1993; Lync 
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and Lorimer, 1994; Taylor et al, 1993; Winkler and Petre, 2007). As with the ID 
simulations, all two-dimensional (2D) simulations studied so far with the state-of-the- 
art neutrino transport (Buras et al, 2003, 2006; Mezzacappa et al, 1998) have yet to 
explode. In their place, 2D gray codes have been used to study how asphericity 
could be generated. Initially the focus was upon neutrino driven convection but 
the anisotropy was found to be small and the pulsar velocities were only of a few 
hundred kms^^ (Burrows et al, 1995; Janka and Mueller, 1994, 1996). While not 
sufficient to explain the observations these studies did note that their simulations 
contained other departures from asphericity: for example, distortions of the stalled 
shock were observed by Burrows et al (1995) who described the shock surface in their 
simulations as botryoidal, i.e. resembling a bunch of grapes. It was further study 
of these distrotions of the stalled shock that have led to the discovery by Blondin 
et al (2003) of a new mechanism for generating asphericity in supernova. Blondin 
et al (2003) found that small aspherical perturbations of a stalled, spherical accretion 
shock could quickly grow and develop large dipole and, to a lesser extent quadrapole, 
modes. This standing accretion shock instability (SASI) will play an important part 
of future supernova simulations and may be a key ingredient in achieving explosions 
based upon first principles. From the multi-dimensional simulations we find that each 
radial slice through the supernova still possess the common features of forward/reverse 
shocks and a global hot bubble region identified from ID but superimposed will be 
a mixture of features typical of asphericity such as additional/internal shocks, local 
bubbles, turbulence, sound waves, etc. Of these additional features in the profile the 
one that has generated most interest is the turbulence. The effects of turbulence in 
supernova profiles was first considered by Sawyer (1990, 1994) and Loreti et al (1995) 
then later expanded upon by Fogli et al (2006a). The correlation function of the 
density fluctuations is an important determinant and a typical first assumption was 
that the fluctuations were 5-correlated. However, it was pointed out by Bcnatti and 
Floreanini (2005) that such a correlation is rather idealised and so recently authors 
such as Fricdland and Gruzinov (2006) and Choubey et al (2007) have switched to 
considering Kolmogorov turbulence instead. Whatever the exact spectrum used the 
general result of all these studies is that turbulence of sufficient strength can lead to 
flavour depolarisation. 

In summary, over the period of just a few years the density profiles used to study 
supernova neutrino oscillations have progressed from simple, static parameterisations 
to profiles taken from aspherical supernova simulations that contain numerous features 
that affect the neutrinos. Later in this review we will examine in greater detail the 
dynamics of the pure MSW effect and how the use of improved density profiles has 
altered our expectations of the information in the neutrino signal. 

1.2. Neutrino self-interaction and collective neutrino flavour transformation 

The contribution from neutrino-neutrino forward scattering, or neutrino self- 
interaction, was noted by some of the first studies of neutrino fiavour transformation 
in stellar collapse (Fuller et al, 1992, 1987) and in the early Universe (Notzold and 
Raffclt, 1988). Shortly after Pantaleone (1992) pointed out the existence of off- 
diagonal elements of the neutrino self-interaction potential in the flavour basis, and 
this observation was then incorporated into the general formalism for neutrino flavour 
evolution laid out by Sigi and Raffclt (1993). From these early studies it quickly 
became obvious that the inclusion of neutrino self-interactions poses a formidable 
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Figure 1. With neutrino self- interaction flavour evolution histories of neutrinos 
propagating along different but intersecting trajectories, e.g. and i/q, are 
coupled because their flavour evolution histories beyond the intersection point, 
Q, will depend on the flavour states of both neutrinos at Q. In addition, the 
flavour evolution histories of non-intersecting neutrino beams such as Vp and Vq 
can also be coupled because the flavour evolution history of Vp beyond P indirectly 
depends on the flavour state of Vq at Q through i/j, . 

challenge (Qian and Fuller, 1995b). The reason for the complexity is illustrated in 
figure 1 . The flavour evolution histories of any two neutrino beams becomes coupled 
due to the self-interaction and, as a result, in order to study the neutrino oscillations 
with neutrino self-interactions one needs to simultaneously follow the flavour evolution 
of all the neutrinos (with different energies, initial flavour states and propagating 
directions) emitted around the same time from all points on the neutrino sphere. We 
must emphasise that coupling of neutrino flavour evolution histories discussed in this 
review is, however, not quantum entanglement. The possible effects of the latter were 
discussed by, e.g. Bell tt al (2003); Fricdland and Lunardini (2003a, b). 

Even for an idealised supernova model with a perfect spherical symmetry one must 
distinguish between the different neutrino trajectories along which neutrinos will have 
travelled different distances for the same radius interval. A simplifying approximation 
(but not necessarily a self-consistent treatment) was proposed by Qian and Fuller 
(1995b) who considered all the neutrino trajectories to be equivalent. In this "single- 
angle approximation" only the flavour evolution of neutrinos along one representative 
trajectory (e.g. the radially directed one) are computed. 

In the two-flavour mixing scheme the flavour state of a neutrino can be represented 
as a classical spin or "magnetic dipole" in the three-dimensional flavour space (Kim 
et al, 1987; Mikheyev and Smirnov, 1986). Using this spin analogy Sigl and Raffelt 
(1993) showed that neutrino self-interaction is represented as coupling between these 
"flavour spins" . (We stress that a flavour spin is a flctitious spin in flavour space and 
that the spin analogy for neutrino oscillations discussed here is not related to flavour 
conversion of neutrinos because of coupling between the neutrino electromagnetic 
moment and a physical magnetic field considered by, e.g. Voloshin et al (1986).) 
While the early studies of neutrinos in supernova (Fuller et al, 1987; Qian and Fuller, 
1995a,b; Qian et al, 1993) treated the self-interaction as an additional matter effect 
in the MSW mechanism, it is not inconceivable that the ensemble of flavour spins 
can exhibit collective behaviors in a similar fashion as spin gases in condensed matter 
physics. In fact, collective neutrino oscillation phenomena have been discussed in 
the context of the early Universe (e.g. Kostclccky and Samuel, 1995). One collective 
oscillation phenomenon which can appear in dense, homogeneous, isotropic neutrino 
gases is that, independent of their energy, all the neutrinos oscillate in the same 
way as a representative neutrino in vacuum (Samuel, 1993). This phenomenon is 
known as the synchronised neutrino oscillation because it can be explained in analogy 
to "a system of magnetic dipoles which are coupled by their self-interactions to 
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form one large magnetic dipole which then precesses coherently in a weak external 
magnetic field" (Pastor and Raffelt, 2002). Another collective oscillation phenomenon 
is known to occur in dense gases with approximately equal number of neutrinos and 
antineutrinos. Both neutrinos and antineutrinos can experience "substantial flavour 
oscillation even for extremely small mixing angles" if the neutrino mass hierarchy 
is inverted (Kostelecky and Samuel, 1993). This behaviour is usually known as the 
bipolar neutrino oscillation because it can explained by using an analogy to a bipolar 
spin system which consists of two coupled and approximately oppositely-oriented 
groups of spins (Duan et al, 2006c). 

Armed with this insight from neutrino behaviour in the early Universe and using 
the single-angle approximation Pastor and Raffelt (2002) and Balantekin and Yiiksel 
(2005) showed that, in a normal neutrino mass hierarchy scenario, the combination of 
synchronisation and the MSW mechanism may indeed cause flavour conversion of both 
neutrinos and antineutrinos if the supernova neutrino luminosities are large enough. 
This phenomenon requires either a large mass-squared difference {Sm? ~ 10 eV^) 
or a dilute supernova envelope so that the synchronised neutrinos encounter a 
MSW resonance near the neutrino sphere. Alternatively, Fuller and Qian (2006) 
envisioned a scenario where the off-diagonal part of the neutrino self-interaction 
potential dominates and both neutrinos and antineutrinos simultaneously achieve 
maximal mixing even for small mass-squared differences. This scenario is essentially a 
special case of synchronisation. Although Fuller and Qian (2006) identified necessary 
conditions for this scenario to occur, it is not clear if a realistic supernova event could 
evolve into this scenario even if these necessary conditions are satisfied. Then Duan 
et al (2006c) showed that a general bipolar neutrino system (e.g. with unequal numbers 
of neutrinos and antineutrinos) can exhibit collective flavour oscillations of either the 
synchronised or bipolar type if the total neutrino flux is larger or smaller than some 
critical value. In the same study Duan et al (2006c) also showed that ordinary matter 
has little effect on collective neutrino oscillations and that bipolar neutrino oscillations 
can occur very close to the neutrino sphere even if the matter density is large. Based 
on their analyses, Duan et al (2006c) mapped out the regimes in supernovae where 
synchronised, bipolar and conventional MSW neutrino oscillations may occur. 

The first multi-angle (i.e. treating different neutrino trajectories as well 
as different neutrino energies) calculations of neutrino flavour transformation in 
supernovae were carried out by Duan et al (2006a, b). The results of their calculations 
show that, contrary to expectations based on the conventional MSW mechanism 
or synchronisation, both neutrinos and antineutrinos can simultaneously experience 
significant flavour oscillations near the neutrino sphere in some inverted neutrino 
mass hierarchy scenarios. Furthermore, when collective neutrino oscillations end, 
neutrinos of different flavours have their energy spectra swapped above a critical 
energy Ef,, a phenomenon now known as "stepwise spectral swap" or "spectral 
split" . With density profiles and neutrino mixing parameters similar to those adopted 
by Balantekin and Yiiksel (2005), Duan et al (2006a, b) observed similar spectral 
swap/split phenomenon with the normal mass hierarchy. They also carried out the 
single-angle calculations with setups similar to the multi-angle ones, and the results 
of single-angle and multi-angle calculations were found to bear many of the same 
qualitative features. This suggests that the single-angle approximation is a valid 
approach to study, at least qualitatively, neutrino oscillations in supernovae. Using 
the single-angle approximation Duan et al (2006b) also showed that the spectral 
swap/split phenomenon can be a result of collective neutrino oscillation mode in which 
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all flavour spins precess collectively. The discoveries made by Duan et al (2006a,b) 
have triggered intense and fruitful studies of collective neutrino oscillations in the 
supernova environment. Similar results were found in calculations using different 
neutrino mixing parameters and matter density profiles (Duan et al, 2007b; Fogli 
et al, 2007) and even for some full three-flavour mixing scenarios (Dasgupta and 
Dighe, 2008; Duan et al, 2008b). Signatures of these oscillation features in supernova 
neutrino signals have been investigated (Chakraborty et al, 2008; Dasgupta et al, 
2008a; Fogh et al, 2008; Gava et al, 2009; Lunardini et al, 2008) and on the theoretical 
side collective neutrino oscillations in a single-angle approximated supernova model 
are now well understood (Dasgupta and Dighe, 2008; Dasgupta et al, 2008c; Duan 
et al, 2007a,c, 2008d; Hannestad et al, 2006; Raffelt and Smirnov, 2007a,b). Collective 
neutrino oscillations in anisotropic environments (Duan et al, 2008c; Esteban-Pretel 
et al, 2007; Hannestad et al, 2006; Raffelt and Sigl, 2007; Sawyer, 2008) and in non- 
spherical geometry (Dasgupta et al, 2008b) have also been examined. Effects of CP- 
violation (Balantekin et al, 2008; Gava and Volpe, 2008) and the case of very large 
matter density (Esteban-Pretel et al, 2008a, b) have also been analysed. The recent 
progress in the understanding of neutrino self-interactions in supernovae has been 
impressive and, like the effect of the dynamic density profiles, has completely changed 
our expectations of the features in the neutrino signal we shall receive from the next 
Galactic supernova. 

1.3. Organisation of the paper 

Due to their importance in the explosion environment the neutrinos and their fiavour 
transformation touch upon many different aspects of supernova phenomenology. For 
lack of space we cannot contemplate a comprehensive review of the entire field so 
we shall focus upon the specific topic of fiavour transformation of active neutrinos 
as they propagate from the PNS to the supernova surface. In doing so we have 
left out many other interesting subjects such as neutrino nucleus interactions (e.g. 
Balantekin and Fuller, 2003), neutrino inelastic scattering (e.g. Martmcz-Pincdo, 2008) 
and the possibility and consequences of sterile neutrinos (e.g. Bcun et al, 2006; Fetter 
et al, 2003; Fuller et al, 2009; Hidaka and Fuller, 2007). We will also not be able 
to address how one can begin to unpick the various signatures we now expect in an 
actual neutrino signal instead referring the reader to the literature (e.g. Jachowicz 
and McLaughlin, 2006; Jachowicz et al, 2008). We begin with some general formalism 
of fiavour transformation of supernova neutrinos in section 2. Historically the effects 
of dynamic matter density profiles and neutrino self-interaction have been studied 
in parallel, i.e. without inclusion of the effects caused by each other. So we will 
independently discuss the MSW flavour transformation with dynamic density profiles 
and collective fiavour transformation in sections 3 and 4, respectively. Finally in 
section 5 we summarise and provide some outlook for where future advances lie. 

2. Neutrino flavour transformation in supernovae: general discussions 

2.1. Neutrino flavour transformation without neutrino self -interaction 

Because neutrinos initially propagate coherently outside the neutrino sphere, the 
fiavour state |'0p('")) of a neutrino with momentum p at point r on its world line 
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can be solved from a Schrodinger-like equation (e.g. Cardall, 2007; Halprin, 1986) 

i^lV'p(r-)) - H{r)\^p{r)) = [i?vac + i?matt(r) + i?..(r)]|7Ap(r)), (1) 

where A is the propagation distance of the neutrino along its world line, H^ac is the 
Hamiltonian for vacuum oscillations, -ffmatt('") is from the forward scattering of the 
neutrino with the ordinary matter, and H^^(r) is the contribution from neutrino self- 
interaction. In equation (1) we have adopted the steady-state approximation. This 
is because neutrinos can traverse regions of interest within very short time, and the 
physical conditions in supernovae can barely change during this period. The flavour 
evolution of a neutrino along its world line can, therefore, be solved from a static 
configuration as outlined in equation (1). We first consider the scenarios that neutrino 
self-interaction can be ignored, i.e. Hiyi,{r) = 0. 

There are two fundamental bases for the neutrinos: the flavour states Iva) 
{a = e, fjL, t) which are the states seen in detectors, and the vacuum mass states 
l^i) (i = 1, 2, 3) which are the eigenstates of i/vac. In the latter basis -ffvac is — up to 
a term proportional to the identity matrix — equal to 

Hi:l^:^diag[mlmlml], (2) 

where the elements of Hitl are defined as {Hill)ij = (i^il^vad^'i), superscript "(v)" 
stands for the vacuum mass basis, E — \p\ is the energy of the neutrino, and rui is 
the mass of Vi. The two bases are related by a unitary transformation 

i 

The unitary matrix U in equation (3) is conventionally parametrised with three 
vacuum mixing angles 612, 023 and ^13, and the three remaining phases 5, ai and 
a2 (Amsler et al, 2008) as: 

U = 



C12C13 S12C13 si3e~"' 

-S12C23 - Ci2S23Sl3e"^ C12C23 - Si2S23Sl3e"' S23C13 
S12S23 - Ci2C23Sl3e"' -C12S23 - 5120233136"^ C23C13 



X diag[e'"i/2^e'"^/2^1], (4) 

where Cij — cos 9ij and Sy — sin 9ij. We will ignore the two Majorana phases ai and 
a2 from this point because they have no effect for flavour transformation of ultra- 
relativistic neutrinos in matter (e.g. Strumia and Vissani, 2006). Various neutrino 
oscillation experiments have measured or put constraints on most of the neutrino 
mixing parameters: Sml2 — rri^ — ra\ ~ Srn^ ~ 8 x 10^^ eV^, |<5r7i23| — |m§ — raW ~ 
5ml^^ ~ 2.4 X 10-3 eV^ sin^ 26li2 ~ 0.86, sin^ 26*23 1 and sin^ 26*13 < 0.19 (Fogh 
et al, 2006b). The CP phase 5 and the sign of 5rn^^ are undetermined. It is 
conventional to call a mixing scheme a normal mass hierarchy if 5rn^^ > and an 
inverted mass hierarchy if ^77123 < 0. 

The Hamiltonian i?matt(''') is a diagonal matrix in the flavour basis 

H^i^^,,{r) = diag[K(r), 0, 0] = ^G^Np^pir) diag[i;(r), 0, 0], (5) 

where superscript "(f)" denotes the flavour basis, Gp — 1.166 x lO^'"^ GeV^^ is Fermi's 
constant, Np^ ~ 6.022 x 10^'^ g"^ is Avogadro's number, p(r) is the matter density, and 
Yeif) is the electron fraction per baryon. We note that the matter temperature outside 
the neutrino sphere is too low to have fi 01 t produced. (However, see, e.g. Dighe and 
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Smirnov, 2000; Esteban-Pretel et al, 2008b for the effects of nonzero "effective" /i- 
abundance in the presence of a very large matter density.) The eigenstates \vf^{r)) 
{i = 1,2,3) of the Hamiltonian H{r) — i/vac + ^finatt(''') constitute the "matter 
basis". We define matter eigenstates \v™'{r)) to be in the same order as that of 
vacuum mass eigenstates \vi), i.e. the eigenvalues fci(r) associated with \v™{r)) satisfy 
fci(r) < k2{r) < fc3(r) for the normal mass hierarchy and A:3(r) < ki{r) < k2{r) for 
the inverted mass hierarchy. Unlike the flavour or vacuum mass states, the eigenstates 
|j'™(r)) as well as the corresponding eigenvalues ki{r) vary with location as the 
electron number density ne(r) = NAp{r)Ye{r) changes. All the above discussions 
also apply to antineutrinos except that there is a sign difference between Hmatti'i') for 
neutrinos and for antineutrinos, i.e. one should take Ve{r) — > —Ve{r) in equation (5) 
for antineutrinos. 

As an alternative to solving equation (1) the evolution of a neutrino can also be 
found by solving for the evolution operator Sp{r,ro). This operator obeys a similar 
equation, 

i-^Sp{r,ro) = H{r)Sp{r,ro). (6) 
and has the additional property that it obeys the product rule: 

Sp{r,ro)^Sp{r,r')Sp{r',r„). (7) 
Equation (6) has a formal solution 

H{r')dX' 



Sp{r,ro) = Texp 



(8) 



where T is the space/time-ordering operator, and Aq, A and A' are the distances along 
the world line of the neutrino that correspond to rp, r and r' , respectively. In practice, 
Spi'f'i fo) is approximated by a sum of truncated series (Kneller and McLaughlin, 2006; 
Kneller and McLaughlin, 2009). 



2.2. The two-flavour approximation 

We have, so far, described everything in terms of three flavours but a quick scan 
through the literature by the reader will reveal that many studies have used only two. 
The reduction in number of flavours comes from observation of how neutrinos evolve in 
the matter basis. The evolution operator Sp in this basis is a diagonal matrix and the 
probability for the neutrino to be in each matter eigenstate is constant except in the 
vicinity of the resonant regions where two of the three eigenvalues ki of H approach 
one another. The evolution of the eigenvalues as a function of the potential is shown in 
figure 2 and from it we see that there are two resonances: the "H" resonance at a larger 
matter density associated with the 6rnl^.^-scale, and the resonance at the JrriQ-scale 
known as the "L" resonance which occurs at a smaller matter density. As a result of the 
large difference between the two measured mass-squared differences, Sthq and Jm^^^, 
neutrino mixing at any resonant location occurs principally between the two matter 
states whose eigenvalues are the closest, and the third matter state remains decoupled 
(Dighe and Smirnov, 2000; Kuo and Pantaleone, 1989b). Therefore, the full three- 
fiavour mixing problem can be approximated as two successive two-flavour mixing 
scenarios that occur at the H and L resonances, respectively. Here the L resonance 
always describes the mixing between |i^™('r)) and \i'2^{r)) but for the H resonance 
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Figure 2. Schematic plots of ki, the eigenvalues of Hamiltonian H, as functions 
of Ve for the normal neutrino mass hierarchy (NH, the left panel) and the inverted 
neutrino mass hierarchy (IH, the right panel), respectively. The "H" and "L" 
MSW resonances are labelled correspondingly. The resonances occur at Ve > 
and Ve < are for neutrinos and antineutrinos, respectively. The eigenvalues of 
the two distinct matter states that are involved in the H resonance arc close to 
each other (but do not equal) at the resonance. 



the mixing states differ being |i^™(T')) and for normal mass hierarchy and the 

states \D™{r)) and |i^™('r)) for inverted mass hierarchy. 

Either the H or L resonance is said to be "adiabatic" or "non-adiabatic" depending 
upon whether the crossing probabiHty Pq is close to zero or closer to unity and 
indeed these are the two natural values. If the electron number density — pY^ 
has a power-law dependence on the distance along the world line of the neutrino, i.e. 
ne(A) oc 1/(A — Aq)", the crossing probability Pq in the two-flavour approximation 
has an analytical solution 

7r7^ F{n, 0^ ) ' 
2 



Pc = exp - 



(9) 



In equation (9) 



= sin^ 29, 



dK(r) 



dA 



is the adiabaticity parameter evaluated at the resonant position r^, where 



5r 



2E 



cos26'v - Ve(r^) = 0. 



(10) 



(11) 



Here 5m? and are the mass-squared difference and effective vacuum mixing angle 
in the two- flavour approximation. The function F(n,9v) in equation (9) depends 
solely on n and 9^ and can be found in Kuo and Pantaleone (1989a) and Kachelriefi 
and Tomas (2001). It is common to find authors using the Landau-Stiickelberg-Zener 
(Landau and Lifshitz, 1977; Stiickelberg, E. C. G., 1932; Zener, 1932) formula for a 
linear density profile for which F = 1. 

When equation (9) is applied to the L resonance with Sm? ~ Sitiq and 9v — 9q 
it turns out that the crossing probability is always very small — as in the sun — and 
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the neutrino transformation is adiabatic (Dighe and Smirnov, 2000). As a result, one 
always has 



Sn 









a* 



H 



and 5*13 



Q!H 

0_ 



/3h 




(12) 



where Sp (Sp) is the matrix of the evolution operator in the matter basis that evolves 
the neutrino from the neutrino sphere to the surface of the supernova, and means 
that the two matrices are different only by some diagonal matrices which do not change 
the probability for the neutrino or antineutrino to be in each matter eigenstate (Kneller 
et al, 2008). In equation (12) an and /3h (Sh and Pn) are the Cayley-Klein parameters 
for S'h (>5h), the evolution operator in the matter basis that describes the two- flavour 
mixing scenario for the H resonance of the neutrino (antineutrino). Up to a common 
phase factor which is irrelevant for neutrino oscillations we can write 



S'h 



an 



and S'h 



an 



Ph 



(13) 



where the dependence of Cayley-Klein parameters on neutrino momentum is implicit. 
The Cayley-Klein parameters satisfy the unitary condition 

|aHp + |/3Hp = l and |«h|' + |/3hP = 1- (14) 

From equations (13) and (14) we see that the crossing probability at the H resonance 
is 



Pi 



H 



1 



(15) 



for neutrinos and 



(16) 



713 



Ph = 1 - I^hI' = I/ShP 
for antineutrinos. 

For the H resonance the appropriate substitutions into equation (11) are 9^ 
and Sm^ ~ zLSml^^, where the plus and minus signs are for normal mass hierarchy 
and inverted mass hierarchy, respectively. But, as figure 2 illustrates, we cannot have 
a situation where an H resonance exists for both neutrinos and antineutrinos. For 
this reason we must have /3h = for a normal hierarchy and /3h = for an inverted 
hierarchy . 



2.3. Neutrino self-interaction potential 

Assuming that the effects of the neutrino medium on the flavour evolution of a test 
neutrino depend on the average flavour content of background neutrinos (Balantekin 
and Pehlivan, 2007; Friedland and Lunardini, 2003a, b), Hamiltonian /f^^ for a 
neutrino with momentum p can be written in the flavour basis as 

Hl'J,p{r) = V2Gf I ^[l-cosdpp,{r)][g,^p,{r) - gl^p,{r)], (17) 

where p and p' are the momenta of the test and background neutrinos, respectively, 
'dpp' is angle between the propagation directions of these neutrinos, and gu{p).p' 
are the matrices of density describing the flavour states of background neutrinos 
(antineutrinos). (In Sigl and Raffelt (1993); Strack and Burrows (2005) the density of 
matrix for an antineutrino is actually defined as [gp{r)]ai3 = [gp.plpa = [Qp.pla/s- This 
definition of the density of matrix with a reversed order in subscripts is helpful if one 
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uses matrices of densities instead of wavefunctions to describe the flavour states of 
neutrinos. See equation (48) and discussions in section 4.2.) Note that we use ^ and 
9 (with or without subscripts or superscripts) to denote angles that are in coordinate 
space and flavour space, respectively. Similarly, the Hamiltonian for an antineutrino 
with momentum p is 

Hl'Ur) = V2Gf J ^[l-cos^pp,{r)][g,^p,ir) - g:^^,ir)]. (18) 

Assuming that the flavour state of each neutrino can be fully described by a 
wavefunction (i.e. without any quantum decoherence) , the elements of the matrices of 
densities for neutrinos and antincutrinos can be written as 

Mr)Uj (19) 
iQ'^,pir)U -E^'^;(^)(^"l^^p(^))(^'^;>WI^/3), (20) 

where v'p (9^) represents a neutrino (antineutrino) beam with momentum p, F^i (^j^i ) (r) 
is its number flux at r, and jV'iy' (p' )) is the flavour state of the neutrino (antineutrino). 
The diagonal elements of a matrix of density give the number densities of neutrinos 
in the corresponding flavours. For example, 

l^[g..p'{r)]ee (21) 

is the total number density of &t r. The off-diagonal elements of a matrix of density 
contain information about neutrino flavour mixing. 

Like in pure MSW flavour transformation, the problem of three-flavour neutrino 
oscillation with neutrino self-interaction is also believed to be factorisable into two 
successive two-flavour mixing scenarios under most circumstances (Dasgupta and 
Dighe, 2008; Dasgupta et al, 2008c; Duan et al, 2008d). In this review will focus on 
the two-flavour mixing scenario at the Srri^^^-scale which is the mostly likely neutrino 
mixing scenario to affect, if at all, supernova dynamics and/or nucleosynthesis. 

3. Neutrino flavour transformation in supernovae: the electron 
dominated regime 

In this section we review the effects of dynamical density proflles on flavour 
transformation of supernova neutrinos. In sections 3.1 and 3.2 we will discuss the 
effects of the forward and reverse shocks and hot bubbles. In doing so we will utilise 
density profiles that are derived from ID supernova simulations and we will focus on 
neutrinos that are emitted radially from the PNS. In section 3.3 we will discuss the 
features in aspherical profiles and flavour transformation of neutrinos that propagate 
along non-radial trajectories. 

3.1. The forward shock 

Until recently any study of supernova neutrino oscillations assumed the density profile 
observed by the outrushing neutrinos to be that of the progenitor star and to be static. 
Roughly this profile obeys an inverse power law i.e. p oc 1/r" with common values of n 
being 2 — 3 (Brown et al, 1982; Dighe and Smirnov, 2000; Kuo and Pantaleone, 1988; 
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Notzold, 1987) though upon closer inspection one finds that profiles of progenitors 
from simulations show deviations from this form (Kneller et al, 2008; Notzold, 1987; 
Takahashi et al, 2003). Whatever progenitor profile used, we must first establish as a 
baseline the effect of the progenitor upon the neutrinos because, a) for the first second 
or so of the neutrino signal this is the profile at the H resonance, and b), it is changes 
from this baseline that reveal the evolution of the density profile. 

An insightful study of the effects of the progenitor was made by Lunardini and 
Smirnov (2003). Using equation (9) with the appropriate mixing parameters for the H 
resonance and adopting an inverse power law profile they distinguished three regimes 
in the parameter space of ^13: 

• The adiabatic regime where 

/ 771 \ 2/3 

sin^ 6*13 > 10-* . (22) 

In this regime 7* ^ 1 and Pa ^ for the normal hierarchy or Pa ^ for the 
inverted hierarchy. 

• The non-adiabatic regime where 

sin2 6'i3 < 10"^ ( — - — ^ ^ . (23) 

In this regime 7^, <C 1 and Pa ^ 1 (or Ph ^ !)• 

• The transition regime where 

In this regime 7* ~ 1 and Ph (or Ph) takes on intermediate values. 

The most significant change in the study of the pure-MSW in supernovae over 
recent years has been the realisation that the profile seen by the neutrinos at later 
times is not the progenitor profile. It is now apparent that prominent features appear 
in the profile of the supernova that are not present in the progenitor and these features 
will affect the neutrinos. The first, generic feature of any supernova density profile is 
the forward shock. Snapshots of density profiles with solely this feature are shown in 
figure 3 which is taken from Kneller et al (2008). The origin of time t is taken to be 
the formation of the shock, an event often called the "bounce" . Observationally the 
bounce would be indicated by the neutronisation burst which occurs a millisecond or 
so afterwards. Kneller et al (2008) mapped into a hydrodynamical code a spherically 
symmetric (ID) density profile resembling the state of the supernova approximately at 
the moment when the forward shock has stalled at r ~ 200 km. The material between 
a 100 km gain radius and the stalled shock was heated in a fashion resembling neutrino 
energy deposition in order to revive the outward motion of the shock. By adjusting the 
energy deposition rate Kneller et al (2008) could control the features of the explosion. 
We should mention that density profiles that only contain a forward shock are not 
typical of iron core-collapse supernova, nevertheless, such a simple profile is a useful 
place to begin to understand the manner in which the explosion affects the neutrinos. 
In the figure the forward shock is the step in density and we can see it racing out 
through the supernova and through the H resonances. For normal shocks the density 
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Figure 3. The density profile of a wealc explosion in a spherically symmetric 
simulation taken from Kneller et al (2008). Displayed are the radial profile at a 
series of snapshot times which are: t = 0.9 s (dotted), 1.8 s (solid), 3.6 s (dashed) 
and 7.2 s (dash dot). The horizontal dashed lines are the H resonance densities 
of, from top to bottom, neutrinos with energies of 5, 10, 20, 40 and 80 MeV. 



(25) 



jump, Apsh/psh, is related to the Mach number, AI, of the shock and the ratio of 
specific heats, Cp/Cy, of the medium via 

Ap,h ^ 2 { AP - 1) 
Psh (Cp/Cy-1) M2 + 2- 
For strong shocks the density jump becomes independent of the Mach number. 

As the shock reaches the H resonance for any given neutrino energy the density 
derivative dp/dr becomes abruptly steeper. This will change the adiabaticity of the 
resonance according to equation (10). If the mixing parameters are such that neutrino 
propagation through the progenitor is adiabatic, i.e. 613 lies in the adiabatic regime 
so that Ph ~ 0, then the arrival of the shock at any given resonance means that 
neutrinos of the same energy now propagate through the profile non-adiabatically, i.e. 
Ph ^1. It is the possibility, first considered by Schirato and Fuller (2002), that Ph 
can evolve with time that allows us to peer inside the supernova as it explodes. Of 
course if 613 is very small so that it lies in the non-adiabatic regime then no change 
will be observed as the shock passes the H resonance region because the H resonance 
is already non-adiabatic for the undisturbed/progenitor profile, i.e. 7* <C 1, Ph ~ 1. 
When the shock arrives all that happens is that 7^ will become even smaller and but 
this has no effect upon Pfj. Such a situation would be disappointing but at least 
we would be able to set an upper limit on ^13 that is competitive with any possible 
terrestrial experiment. So for the rest of this discussion of the dynamic MSW effect we 
shall consider only the case where ^13 lies in the adiabatic regime, i.e. 6*13 is such that 
equation (22) is satisfied, and more specifically adopt a value sin^ 26*13 = 4 x 10~*. 

With our understanding of how Ph can change as the star explodes we can begin 
to map out our expectations for the evolution of Ph as a function of time and energy. 
Let us consider the particular example shown in figure 3. For the first couple of 
seconds the profile at the H resonances of neutrinos with energy 5 — 80 MeV is still 



Neutrino flavour transformation in supernovae 16 




E [MeV] 

Figure 4. The crossing probability as a function of neutrino energy for a 
supernova simulation containing only a forward shock. Density profiles from this 
simulation are shown in figure 3 and the snapshot times are: t = 3.0 s (solid), 
3.5 s (dashed), 4.0 s (dash-dot) and 4.5 s (dash-double dot). The figure is taken 
from Kneller et al (2008). 




Figure 5. The crossing probability as a function of time for selected neutrino 
energies as they pass through the density profiles taken from the weak explosion 
shown in figure 3. The energies are: i? = 10 MeV (solid), 15 MeV (dashed), 
20 MeV (dash dot), 25 MeV (dotted) and 30 MeV (dash double dot). The figure 
is taken from Kneller et al (2008) 



the undisturbed progenitor and only after this initial delay does the shock appear in 
the H resonance region. Then si t — 2.0 s the shock will change the adiabaticity of 
the 5 MeV neutrinos but will not begin to affect the 40 MeV neutrinos until later at 
3.6 s because it takes longer for the shock to reach their H resonance density. Also 
around 3.5 s the shock ceases to intersect the 5 MeV resonance density so for this 
energy the adiabaticity will return to Ph ~ 0. Only later at t = 7.2 s does the shock 
cease affecting the 40 MeV neutrinos. 

In figures 4 and 5 we show the actual numerical results for the crossing probability 
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Ph as a function of energy at various snapshot times, and as a function of time for 
fixed neutrino energies respectively as neutrinos pass through density profiles taken 
from the weak supernova shown in figure 3. We clearly see the arrival of the shock at 
early times as the sudden jump in Pjj at low energies which then sweeps up through 
the spectrum. Note how the width of the shock feature increases with time. We also 
see that higher energy neutrinos are affected after the lower energies and that the 
effect is of longer duration as expected. From figures 4 and 5 we can identify three 
basic quantities that describe the evolution of 

• Esh{t), the highest energy affected by the shock at any given moment t, 

• AEshit), the width of the shock feature in the spectrum at any given t, and 

• Atsh{E), the time period over which the shock affects neutrinos of energy E. 

From figures 4 and 5 we see that both -E'sh(i) and A_Esh(^) increase with t and that 
Atsh{E) increases with E. This behaviour is easily understood. 

The highest energy affected by the shock, Esh{t), at time t corresponds to the 
moment when the bottom of the shock reaches the resonance density for that particular 
energy as shown in figure 3. Esh{t) is a monotically increasing function of t which is 
a generic prediction and occurs simply because the resonance density is inversely 
proportional to the neutrino energy, the density profile is a monotonically decreasing 
function of distance and because the shock wave is generated at the core of the star 
where the density is highest. 

The width of the shock feature, AEsh{t), also monotonically increase with t. This 
can be understood as follows. The resonance condition (11) allows us to relate Esh{t) 
to the value of the neutrino potential at the bottom of the shock, Vsh{t), by 

_ Sm^ cos 20V 

^sh - TpF^ • (2bj 

The lowest energy affected by the shock, E^h ~ AEsh are those neutrinos whose 
resonance density is the top of the shock, i.e. Vgh + AVsh where AT4h is the jump 
in the neutrino potential across the shock. Using equation (26) and exploiting the 
fact that AVsh/V^h = Ap^h/psh — where Pshit) is the density at the bottom of the 
shock and Apsh(t) the density jump — we can derive the simple relationship that 

APsh ^ Apsh ^27) 

Esh — AEsh Psh 

Note that this is independent of the mixing parameters. If the ratio of specific heats 
of the material does not change rapidly with distance then Apsh/Psh is constant for 
a strong shock — equation (25) — so the ratio AEsh/{Esh — AEsh) is also invariant. 
This relationship shows us that AEsh oc Esh- Thus, if Esh monotonically increases 
with t then so must AEsh- 

Finally, that Atsh^E) increases with E is also just a consequence of Esh increasing 
with t- The density psh{t) is the point where the shock attaches to the progenitor 
profile so if we assume a form for the density profile of the progenitor then we can 
relate Psh(^) to the position of the shock, rghit). For illustrative purposes we shall adopt 
for the profile the power law p — Ci, (r^/r)" with C^, the constant of proportionality 
and some scale but it is possible to do better if the progenitor star can be identified. 
Using the adopted profile we find that Esh is related to the shock radius by 

fVSGFY.C.EshV^" 

rsh = T-^ ^ ■ 28 

\ niffdm'^ cos26'v / 
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At a time t' the shock will have moved outwards through the star and will be located 
at a position r'^^ = rsh(i'). As it does so the shock feature will have swept up 
through the neutrino spectrum with E'^^ = Esh(t') now the highest affected energy 
and the spectral width will now be AiJ^j^ — AEsh{t'). From equation (28) we see that 
fsh/r'g^ = {Ye Esh/Y^ ^'shY^^ where is the electron fraction at r^j^. If we select t' to 
be the moment when neutrinos of energy Egh cease to be affected by the shock then 
we must have Egh = E'^^^ — AE'^^^. The time difference t' — t = At{Esh)- If Vsh is the 
average shock velocity over this time period then we can derive that 

1/" /vf\ 1/" 



771/ \ i/" /Vl 



1 



(29) 



The term in square brackets is a constant for a constant Apsh/Psh and a medium with 
constant thus At{E^\^ is proportional to r^j^. With this discovery we immediately 
see that if r^h increases with Esh_, so must At{Esh) increase with Esh_. 

While equation (29) shows us that At{Esh) must increase with i?sh perhaps what 
is more interesting is that we can turn the equation around and use the observations 
of Esh{t), AEsh{t), and Atsh{E) to determine the average velocity of the shock as it 
travelled from to r'^^; that is 



?'sh(£'sh) 

Vsh - 



Atsh{Esh) 



1/n /vf\ 1/" 

- 1 



(30) 



As an alternative, we could simply exploit equation (28) to relate the observed Esh{t) 
to construct rghit) which can then be fit with a parametric from. The shock velocity 
would then be found as the derivative of rsh(i) but there is additional information in 
the fit. The shock is formed at the PNS and current models of supernova indicate 
it stalls at r ~ 200 km. Somehow the outward motion of the shock is revived and 
thereafter it races through the mantel of the supernova. The time spent by the shock 
in the stalled position is not known but suggestions typically lie in the range of ^ 0.5 s. 
The velocity of the shock through the star is also not known but a reasonable initial 
guess is that the shock velocity is approximately constant because the profile is close 
to a form (Sedov, 1959). Thus our expectation might be that the position of the 
shock after its revival is a simple linear model given by 

rsh{t) = rstall + Vsh {t - tstall) (31) 

where rgtaii is the radius at which the shock stalls and ^staii is the time at which the 
shock was revived. The shock velocity v^h is of course the gradient of equation (31) 
and the intercept is equal to rgtaii — ^sh ^staii- So if we set rgtaii to zero then we derive 
that tstaii — '"sh/wsh and obtain the earliest time that the shock can have been revived. 
Thus, if the neutrino signal from the next Galactic supernova is of sufficient quality we 
might be able to determine a lower bound for tstaii and in doing so we would be able to 
confirm or refute a fundamental component of the core-collapse supernova paradigm: 
that the shock stalls and is revived. 



3.2. Hot bubbles and reverse shocks 

The density profiles of the supernova simulation used in the previous section were from 
a very weak supernova simulation and contained only a forward shock. The shock was 
revived by heating the material above the PNS but the energy deposition does not not 
switch off once the forward shock has been revived but rather will continue for some 




Figure 6. Density profiles containing a reverse sliock. Thiis figure is adapted 
from a figure in Kncllcr et al (2008). Tiie snapshot times are: t = 1.5 s (daslied), 
2.0 s (dot double-daslied), 2.5 s (dot dasiied) and 3.0 s (solid). 



time after. The continued heating creates a wind that pushes outwards to create a 
hot bubble behind the forward shock. It was actually profiles of this type that were 
studied by Schirato and Fuller (2002). If we up the rate of energy deposition then 
strength of the wind can grow to the point where its velocity becomes greater than the 
sound speed. This situation leads to the formation of a second shock which faces the 
PNS rather than the exterior (Burrows et al, 1995; Janka and Mueller, 1995), i.e. it 
is reversed compared to the forward shock. Profiles with reverse shocks are shown in 
figure 6. The behaviour of the reverse shock has been studied by Tomas et al (2004), 
Arcones et al (2007) and Kneller et al (2008) in ID and 2D using hydro codes of 
varying degrees of sophistication. It has been found that its behaviour is much richer 
than the forward shock due to its greater sensitivity to the mechanism by which it 
was formed, i.e. the energy deposition. If the heating is sustained for a long period 
then the reverse shock trails the forwards shock through the star, but if the energy 
deposition is over a briefer period then the reverse shock has been seen to diminish 
in size, stall and then to turn around and head back to the core. This can be seen in 
figure 6. Also we see from the figure how the densities either side of the reverse shock 
will change in relation to the forward shock as a function of time. 

The hot bubble and the reverse shock also affect the neutrinos. If a profile contains 
either feature then we see that, at any given moment, there are some neutrino energies 
that will experience three (or more) resonances. In order to compute the effects of 
multiple resonance one must be careful less certain effects are ignored. Examples 
of calculations of Ph as a function of neutrino energy through supernova profiles 
containing reverse shocks are shown in figures 7 and 8. While some parts of figure 7 in 
particular resemble figure 4 we see also that something has noticeably changed. The 
crossing probability is seen to oscillate rapidly for some ranges of neutrino energies 
and the reason is due to "phase effects". That phase effects are present in Ph for 
supernova profiles was first noticed by Fogli et al (2003). They were found again by 
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E [MeV] 

Figure 7. The crossing probability as a function of neutrino energy for 
a supernova simulation containing a reverse shock. Density profiles for this 
simulation are shown in figure 6. The snapshot time is t = 2.0 s and the 
energy resolution is 100 keV. The figure is adapted from KncUer et al (2008) and 
corresponds to the ID model with a total energy deposition Q = 3.36 X 10^^ erg. 




E [MeV] 

Figure 8. The crossing probability as a function of neutrino energy for a 
supernova simulation containing a reverse shock. Density profiles from this 
simulation are shown in KncUer et al (2008) and correspond to the model where 
the total energy deposition is Q = 4.51 X 10^^ erg. The snapshot time is f = 2.0 s 
and the energy resolution is 100 keV. 

Knellcr and McLaughlin (2006) who also presented some basic arguments for why the 
phase effects appear. The phase effects and their detectability were then the focus of 
Dasgupta and Dighc (2007). 

While, perhaps, initially surprising the origin of the phase effects is quite 
simple. The neutrino wavefunction after passing through multiple H resonances is 
'0(r^) = S-a_{ri,, Ru)4'{Rv) where r^, is some point between the H resonances and the 
L resonance, and Ri, is the radius of the neutrino sphere. We can factorise S^i by 
dividing up the profile into regions within which there is just one H resonance, i.e. 
just one location in each region where equation (11) is satisfied. Using the product 
rule, equation (7), Sy{ is the product of evolution operators over the sub-domains: 
Sy{ = ... S'h2 'S'hi where S-a_i is the S matrix for passage through the i'th H resonance 
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encountered by the neutrino. All these Sm matrices, up to a phase which is irrelevant 
for neutrino oscillations, have the same structure in the two flavour approximation, 



(32) 



where am and f3m are, again, Cayley-Klein parameters. The crossing probability 
for the i'th resonance is Pm — l/^HiP = 1 ^ |aHjP- Even for the case of three 
H resonances we obtain a long and complicated expression for the net effect of the 
multiple resonances so to keep things simple we shall consider the example of just two 
resonances. In this case the net H resonance crossing probability, Ph, is just 

Ph = Pm (1 - Pm) + (1 - Pm) Pm + 2^{a^iaml3mPm) (33) 

= P2 (1 - Pi) + (1 - P2) Pi + 2VP1P2 (1 - Pi) (1 - P2) COS0H (34) 

where (/)h is a phase formed from the phases of the a's and /3's. The first two terms 
are what we would expect if there were no correlations; the last term represents 
interference and depends upon the relative phases of the a's and /3's, hence the name 
"phase effects". Note the amplitude of the oscillatory term is always smaller than 
the constant and maximal when Phi = Ph2 = 1/2. From this result we can begin 
to understand the results of Tomas et al (2004) for profiles containing both reverse 
and forward shocks. The net effect upon the neutrinos of a given energy depends 
upon the relationship between the shocks and the resonance density. If the two shocks 
intersect completely different resonance densities then we have a situation where, at 
most, either \/3i\ — 1 or |/32| — 1 but not both. But if we have a case where the 
two shocks affect the same resonance densities then we can have a situation where 

— 1^2] — ^ and the two non-adiabatic transitions cancel. As we see in figure 6, the 
two shocks do not always overlap completely and so, due to the relative proportions 
of the forward and reverse shocks, the cancellation may be over just a limited range 
of neutrino energies while energies slightly larger or smaller may experience just one 
shock. This can create what Tomas et al (2004) referred to as the "double dip" in the 
crossing probability Ph. 

But in the more general case with many H resonances some of which are neither 
exactly adiabatic nor non-adiabatic we must consider the effect of the oscillatory terms 
such as the one in equation (34). Dasgupta and Dighe (2007) show that the phase (p^ 
appearing in equation (34) is approximately 

where the rmiE) are the (energy dependent) positions of the resonances i.e. those 
locations that satisfy equation (11). Equation (35) is just the integral of the difference 
between the eigenvalues over the distance between the two resonances. In order to 
compute it properly one needs the density profile Ve{r) but for our purposes we 
approximate the result by neglecting the first term under the square root since we are 
in the vicinity of the resonances where this term vanishes. Thus 0h is approximately 

|(Sm2| sin26'vPH 

'^^ " 2E 

where Lh = '"H2 — ''hi is the distance between the resonances. Typically this distance 
is much greater than the wavelength of the oscillations (~ E /6m^ sin29v), so that (/)h 
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is a large number (0h ^ 1). Now let us consider a change in E and/or Lh- The 
resulting phase change (50h is 

"(5Lh SE' 



Lh L; 



(37) 



We require just a change (5(/)h — 2t: in order for the oscillatory term in equation (34) 
to cycle through one period. If 0h is large, say (/)h ~ 10'^, then a fractional change in 
the energy of just SE/E ~ 10~^ is enough to give us J0h ~ 27r. For neutrinos with 
energy E ^ 10 MeV we see that Ph will cycle over one period with a change of energy 
of just 5E 10 keY. 

From our understanding of equation (34) we can begin to digest the results for 
Ph as a function of neutrino energy shown in figure 7 and 8. In the first figure 
we again observe the presence of the forward shock in the spectrum because the 
crossing probability makes the transition from adiabatic to non-adiabatic propagation 
at 50 MeV but in this same figure we also see that below E ^ 20 MeV the 

crossing probability starts to oscillate as a function of E: this occurs because those 
energies now passes through multiple H resonances and so phase effects appear. In 
figure 8 phase effects dominate the entire spectrum. In both cases the oscillations are 
very rapid which is simply an indication that the resonances are spaced by a large 
number of oscillation wavelengths. 



3.3. Aspherical and turbulent profiles 

As we indicated back in section 1 there is now ample evidence that supernovae are 
aspherical and there has been significant recent progress in identifying what processes 
may be the root of the asphericity. The standing accretion shock instability was 
found by Blondin et al (2003) to generate large dipolar and quadrapolar modes from 
small perturbations of a stalled, spherical accretion shock. When the simulations 
were repeated in three-dimensional (3D) differential, post-shock fiow was found 
indicating the "spinning up" the PNS (Blondin, 2005). More recent work (Blondin and 
Mezzacappa, 2006, 2007; Blondin and Shaw, 2007; Iwakami et al, 2008; Ohnishi et al, 
2008; Scheck et al, 2006) have confirmed the result and furthered the understanding 
of the SASI. 

As an example of an aspherical density profile we show in figure 9 a snapshot, 
taken from Kneller et al (2008), of an aspherical explosion. For this calculation the 
same "initial" density profile used in ID was mapped into a 2D hydrodynamical code 
and heated to revive the motion of the forward shock. The total energy deposition 
was Q « 3 X 10^^ erg, the canonical explosion energy for the matter portion of 
the supernova. To generate the asphericity the heating was inhomogeneous with 
preferential energy deposition along the equator. From the figure we see that the 
density profile along each line of sight can be very different but within the density 
profile along any given radial slice we can still identify the forward and reverse shocks 
and a hot bubble region between them. The radial positions of the two shocks now 
depends upon the polar angle i9. Upon closer inspection we find even in 2D that the 
reverse shock also abates and moves back to the core. Superimposed upon the generic 
template are "fluctuations" generated by more shocks, localised bubbles, sound waves, 
etc. 

Snapshots of Ph as a function of neutrino energy are shown in figure 10. The 
transition from the progenitor profile to the forward shock appears once again as the 
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Figure 9. The density as a function of the radius and angle in a 2D supernova 
model at t = 2.5 s taken from KncUer et al (2008). The forward shock is located 
to the left of "A" , the reverse shock is the step-up in density found to the right of 
"B" , and one of the many localised features in the profile is to the right of "C" . 

transition to Pa ~ 1 and we also notice how Ph drops suddenly around E ^ 10 MeV 
for t = 2.4 s when the profile becomes adiabatic again. Below E' ^ 10 MeV 
phase effects appear as the profile develops multiple resonances for these neutrino 
energies. Given the complexity of these profiles and our understanding of how multiple 
resonances affect neutrinos it comes as little surprise then that phase effects are 
prominent. By t — 5.4 s they dominate the entire spectrum above 20 MeV. After 
t = 6.4 s the forward shock, hot bubble, etc. have largely swept through the spectrum 
and we are able to notice a narrow range of energies surrounding ~ 15 MeV where 
phases effects again appear moving down through the spectrum. This is the signature 
of the reverse shock returning to the core. 

The crossing probabilities shown in figure 10 are very similar to those shown 
earlier in figures 7 and 8. After comparing figure 10 with those figures one has a sense 
that the phase effects are somewhat stronger in 2D. In fact the phase effects are so 
strong that the crossing probability over wide swathes of the neutrino energy range 
appears to be essentially random with, more-or-less, uniform distribution from zero 
to unity. This impression of the uniformity of the distribution of Pjj is confirmed by 
figure 11 where we plot the histogram of Pr from the middle panel of figure 10 for 
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Figure 10. The H resonance crossing probability Pjj a-s a function of neutrino 
energy for a radial slice at i? = 25° through the 2D supernova model. From top 
to bottom the snapshot times are t = 2.4 s, 5.4 s and 6.4 s. The figure is taken 
from Kneller et al (2008). 

neutrino energies between 30 MeV and 80 MeV. While not proof, figure 11 hints that 
the density profiles from 2D supernova simulations may cross over into the regime of 
turbulence. The effects upon neutrinos of a turbulent, noisy, density profile either 
in general or upon the solar density profile date back to Schafer and Koonin (1987), 
Krastev and Smirnov (1989), Sawyer (1990), Loreti and Balantekin (1994), Balantekin 
et al (1996) and many others thereafter. Turbulence in supernova was first considered 
by Loreti et al (1995) and then later by Fogli et al (2006a) in light of the fact that the 
evolving density profile was found to leave its imprint on the neutrino signal. In both 
cases (5-correlated density fiuctuations were assumed. More recently Fricdland and 
Gruzinov (2006) and Choubcy et al (2007) have used Kolmogorov fiuctuation spectra. 
Whatever the fiuctuation spectra used the same result emerges: if the turbulence is 



Neutrino flavour transformation in supernovae 



25 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



"0 0.2 0.4 „ 0.6 0.8 1 

Ph 

Figure 11. A histogram of Ph for neutrino energies between 30 and 80 MeV 
for neutrinos passing through the t = 5.4 density profile shown in figure 10. The 
histogram is formed from 501 points and there are 20 bins so a uniform distribution 
would give an average count of 25 in each bin. 

strong then Ph becomes a random variate uniformly distributed over the range [0, 1]. 
While this prediction seems to be in general agreement with figure 11 we remind the 
reader that what is plotted is the histogram of Ph for different neutrino energies at 
one snapshot time which is something different. For Fogli et al (2006a) fluctuations 
were restricted to a scale of less than 10 km while the spectra used by Friedland 
and Gruzinov (2006) permitted fluctuations on much larger scales. Friedland and 
Gruzinov (2006) showed that it is the density fluctuations on the scale of the oscillation 
wavelength at the resonances that contribute most and that the amplitude of the 
density fluctuations need only be greater than 5p/p > O.l(0i3)^/'^ in order to enter 
the strong turbulence limit. For the value of ^13 used when computing the results in 
figure 10 we require fluctuations of just 1% to reach this limit so indeed turbulence 
should be expected. Detecting turbulence in a neutrino signature would be a clear 
signature that the supernova exploded aspherically but one must not forget that we 
shall receive this signal along just one line of sight which makes demonstrating that 
the signal is the result of turbulence quite difficult. 

4. Neutrino flavour transformation in supernovae: efltects of neutrino 
self-interaction 

In this section we review neutrino oscillations with neutrino self-interaction in 
supernovae. For lack of space, we will focus on the key results of the two-flavour 
numerical simulations of supernova neutrino oscillations by Duan et al (2006a,b) 
and the analytical explanations of the results obtained under the single-angle 
approximation. In section 4.1 we describe in detail the two major approximations, 
i.e. the single-angle and multi-angle approximations, that are now employed to treat 
neutrino oscillations with neutrino self-interaction in supernovae. We also highlight 
the key results in the numerical simulations carried out by Duan et al (2006a, b). The 
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Schrodinger-like equation that governs neutrino flavour transformation, although well 
suited for studying flavour evolution of a single neutrino, is inconvenient when the 
neutrino potential Hi,^ is not negligible. In section 4.2 we introduce the "language", 
the notation of neutrino flavour isospin (NFIS), that we will use to analyse and 
visualise collective neutrino oscillations. In section 4.3 we describe an adiabatic MSW- 
like flavour evolution of dense neutrino gases in the presence of ordinary matter. This 
type of neutrino oscillation echos the early explorations (e.g. Qian and Fuller, 1995b) 
where neutrino self-interaction was included as an additional matter effect in the 
conventional MSW mechanism. In section 4.4 we analyse the stability of bipolar 
neutrino systems such as supernova neutrinos by utilising the pendulum analogy as 
well as the simple concept of energy conservation. In section 4.5 we consider the 
matter effects and try to gain some insights into the qualitative behaviour of neutrino 
systems as neutrino number densities decrease. We explain why supernova neutrinos 
do not follow the MSW-like flavour evolution all the way through. In section 4.6 we 
discuss a collective neutrino oscillation mode which is believed to cause the spectral 
swap/split, a novel phenomenon revealed in numerical calculations. In section 4.7 we 
briefly report current understandings of collective neutrino oscillations with full three 
flavours and/or in realistic supernova environments which are highly inhomogeneous 
and anisotropic. 

4-.1. Approximations and numerical results 

It is clear that neutrino oscillations in the collective regime must be treated in a 
way different from that in the pure MSW regime. In the pure MSW regime, the 
flavour evolution of any single neutrino can be calculated without knowing the flavour 
states of other neutrinos, as long as the matter densities along the world line (i.e. the 
matter profile) of this neutrino is given. For a complicated 3D matter profile, the 
flavour evolution histories of neutrinos propagating along different trajectories can be 
different even if they are initially in the same flavour state and have the same energy. 
In the pure MSW regime, this means that the same algorithm for calculating flavour 
evolution of a single neutrino needs to be run repeatedly for each distinct neutrino 
beam. In the collective regime, however, because the flavour evolution histories of 
neutrinos propagating along different trajectories are all coupled, the flavour states 
of all distinct neutrino beams must be followed simultaneously. While the former 
problem is ready to be solved given enough computing time, the latter poses such a 
great numerical challenge that it has not been tackled yet. 

Partly because of the complexity of the problem, the supernova models adopted so 
far to treat neutrino flavour transformation with neutrino self-interaction are all simple 
and ideal. In these models the supernova environment is spherically symmetric and 
neutrinos and antineutrinos are emitted isotropically from an inflnitely thin neutrino 
sphere with radius R^. Neutrinos and antineutrinos are in pure flavour states at the 
neutrino sphere and encounter only forward-scattering with other particles outside 
the neutrino sphere. The relativistic effects such as redshift of neutrino energies and 
gravitational bending of neutrino trajectories are ignored. 

In such supernova models nonequivalent neutrino trajectories at radius r can 
be distinguished by -d, the angle between the radial direction and the propagation 
direction of the neutrino. Equation (17), therefore, becomes 
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/ dE'f,^,iE'){Dp\i;,^,{r)){^P,^,{r)\D^) , (38) 

where 'i?max(?') = arcsin(i?y/r) is the maximum possible value for i9 at r, rt*°*(r) is the 
total number density of neutrinos and antineutrinos at r, ,) and /^^(p ,)(-£') are 

the fractions of total number flux and normalised energy distributions of neutrinos 
(antineutrinos) with flavour a' at the neutrino sphere, and IV'i/^, (!/„,)('')) is the flavour 
state of a neutrino (antineutrino) initially in flavour a' and (implicitly) with energy 
E' and trajectory angle -d' . 

Even with great simplifications, the "multi- angle approximation" or treatment 
discussed above is hard to analyse and can be computationally intensive. For 
these reasons, the early works on the subject adopt what is called the "single- 
angle approximation" . Under the single-angle approximation neutrino trajectories 
with different i9 are assumed to be equivalent, and the flavour evolution of neutrinos 
propagating along a representative trajectory (e.g. the radial trajectory) are computed. 
The single-angle approximation further simplifies equation (38) and one obtains 

{Hi'Xr)U^^\j2^^^, r AE'UjE')M,^,{r)){^,^,{r)\vp), 

/ AE'f,^,{E'){vp\4>,^,{r)){i^,^,{r)\v^) , (39) 



where 



^i{r) = 2\/2GF<°'(r)C(r) (40) 

is the "effective strength" of neutrino self-interaction at r. In equation (40) C{t) is a 
geometric factor that takes partly into account the angle effects and is 

/1._m(1-cosWos^')_1 



C{r) 



1-1 1- 



(41) 



, , ,d(cosi?') 2 

if the representative neutrino trajectory is radially oriented. (See Dasgupta et al, 
2008b for a discussion of the single-angle approximation in non-spherical geometry.) 

Figure 12 shows the results of the first numerical simulations by Duan et al 
(2006a, b) of flavour transformation of supernova neutrinos employing the multi-angle 
treatment. (Also see Duan et al, 2008a for the movies for the entire simulations.) For 
comparison the results of the corresponding single-angle calculations are also shown 
in the same flgure. These calculations adopted two-flavour mixing schemes with mass- 
squared differences Srn^ ~ ztSm'^^^ and a small mixing angle (^v = 0.1) representing 
the unknown value of ^13. The mixing occurs between the electronic flavour and the 
r' flavour, a linear combination of the muon and tau flavours. The matter proflle and 
the neutrino fluxes/spectra used by these calculations correspond to a late-time epoch 
when r-process nucleosynthesis is supposed to occur. 

Figure 12 shows that, in both the normal and inverted neutrino mass hiearchy 
cases, neutrinos and antineutrinos with various energies propagating along different 
trajectories can simultaneously experience signiflcant flavour transformation near 
the PNS. This is clearly different from a standard MSW-flavour-transformation 
phenomenon. Interestingly, although flavour evolution histories of neutrinos 
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Figure 12. Energy- averaged survival probabilities {P^u) for (left panels) 
and Pe (right panels) as functions of radius r for the normal (upper panels) and 
inverted (lower panels) neutrino mass hierarchies, respectively, when neutrino 
self-interaction is taken into account. The solid and dashed lines give average 
survival probabilities along the radial and tangential trajectories, respectively, 
as computed in the multi-angle simulations. The dotted lines give the average 
survival probabilities computed in the single-angle simulations. Figure adapted 
from Duan et al (2006a). 



propagating along various trajectories are different in multi-angle calculations, the 
single-angle calculations seem to share some key qualitative features with their multi- 
angle counterparts. 

Figure 13 reveals an interesting result of collective oscillations of supernova 
neutrinos which is obtained using the multi-angle treatment. It shows that, when 
neutrino fluxes have dropped off and collective oscillations end, survival probability 
Pun for the neutrino becomes approximately a step function of neutrino energy E^. 
In addition, the directions of these step functions are opposite for the normal and 
inverted neutrino mass hierarchy cases. The corresponding single-angle calculations 
yield results that are qualitatively consistent with those employing the multi-angle 
treatment. The phenomenon seen in figure 13 is called "stepwise spectral swapping" 
because v^^s and Vr'^s apparently swap energy spectra at energies below (above) a 
critical energy in a normal (inverted) neutrino mass hierarchy case (Duan et al, 
2006b). This phenomenon is also known as the "spectral split" since E^, "splits the 
transformed spectrum sharply into parts of almost pure but different flavours" (Raffelt 
and Smirnov, 2007b). 

Neutrino flavour transformation in supernova with neutrino self-interaction has 
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Figure 13. Results from the same multi-angle simulations as shown in figure 
12. This figure shows survival probabilities P^^ for neutrinos (left panels) and 
antineutrinos (right panels) at radius r = 225 km as functions of both neutrino 
energy Ei, and emission angle 'Sn^ , the angle between the propagation direction 
of the neutrino and the normal to the neutrino sphere. The upper panels employ 
a normal neutrino mass hierarchy, and the lower panels employ an inverted 
neutrino mass hierarchy. The colour scale in the plot denotes the neutrino survival 
probability with P^^ = 1 being the lightest. Figure adapted from Duan et al 
(2006a). 

since been studied using both the single-angle and multi-angle treatments. All these 
calculations confirm that single-angle and multi-angle calculations do share the key 
common features shown here. Further study by Duan et al (2007b) showed that the 
swap/split energy Ef, decreases significantly with 6^ in the normal mass hierarchy case 
and is not sensitive to 9^ in the inverted mass hierarchy case. Adopting a different 
matter profile with a thick supernova envelope, Fogii et al (2007) found no significant 
flavour transformation of supernova neutrinos in the normal mass hierarchy case and 
a similar spectral-swap/split phenomenon in the inverted mass hierarchy case. 

The single-angle and multi-angle treatments are clearly not equivalent. In fact, 
the single-angle approximation is not even a self-consistent treatment of the problem. 
For example, neutrinos propagating along trajectories with different values of travel 
different distances for any given radius interval. Therefore, if one uses the results from 
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a single-angle computation and calculates flavour evolution for neutrinos propagating 
along a trajectory that is different from the representative one, one expects to find 
results that will contradict the single-angle assumption. 

Meanwhile, it is also clear that the single-angle treatment is a much simpler 
model than the multi-angle one. Although the multi-angle treatment approaches 
the inhomogeneous and anisotropic nature of the supernova environment in a self- 
consistent way, it is very difficult to study analytically. Under the single-angle 
assumption, however, supernova neutrinos are essentially treated as a homogeneous, 
isotropic neutrino gas with effective neutrino number density 

nfir) - <°*(r)C(r) (42) 

that expands with "time" r. Because the results of the single-angle and multi- 
angle numerical calculations available so far share many qualitative features, Duan 
et al (2006b) conjectured that collective neutrino oscillations in supernovae can 
be understood, at least qualitatively, by studying similar phenomena in isotropic, 
homogeneous neutrino gases. 

For this reason we shall discuss collective flavour transformation of supernova 
neutrinos under the single-angle treatment or, equivalently, a homogeneous, isotropic 
neutrino gas which expands with "time" r. In the following discussions neutrino 
number densities must be understood as the effective ones with the geometric factor 
C(r) included. [For truly homogeneous, isotropic neutrino gases one has C(r) = 1 and, 
therefore, n'^^{r) — 7i*°*(r).] These analyses shall offer valuable insights into the full- 
fledged problem of collective neutrino oscillations in realistic supernova environments. 
In section 4.7 we will briefly summarise the current understandings of collective 
neutrino flavour transformation in inhomogeneous, anisotropic environments. 



4-. 2. Neutrino flavour isospin 

Because Hi,^, is essentially a sum of matrices of densities for all background neutrinos 
and antineutrinos, it is convenient to use matrices of densities instead of wavef unctions 
to describe the flavour states of neutrinos when neutrino self-interaction is important. 
For a two-flavour mixing scheme, a matrix of density 

g,y,E ^ ^[n„,E + ^ ■ Pi^,e] (43) 

is equivalent to polarisation vector P^^e because the trace of g^^E is not relevant for 
neutrino oscillations. In equation (43) we have used energy E to identify a neutrino (or 
antineutrino) mode in a homogeneous, isotropic neutrino gas, and n^^E is the number 
density of neutrinos with energy E. Note that we use "X" to indicate a vector in 
flavour space (as compared to vector "X" in coordinate space), and a is such a vector 
consisting of the three Pauli matrices. Viewing polarisation vectors as "magnetic 
spins" and neutrino self-interaction as coupling between these spins. Pastor et al 
(2002) elucidate the physics mechanism behind synchronisation, a collective neutrino 
oscillation phenomenon discovered in earlier numerical simulations (Samuel, 1993). To 
fully exploit this spin analogy, we will adopt the notation of neutrino flavour isospin 
(NFIS) deflned by Duan et al (2006c). 

For a neutrino (antineutrino) state described by wavefunction '<piy(i^),Ej the 
corresponding NFIS is defined as 

^ I V'I.bIV'.^.b for neutrino, ^^^^ 

I {o'y4'v,E)'^ ^{cryipi^.E) for antineutrino, 
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where 

with the plus and minus signs for the neutrino and the antineutrino, respectively. 
The cTy-transformation of ipp,E in equation (44) eliminates the superficial distinction 
between neutrinos and antineutrinos in the two-flavour mixing scheme. In this 
notation flavour states of neutrinos and antineutrinos are represented by spins in 
flavour space with different values of oj. In the flavour basis (i.e. ip^^E = {i'a\i^u,E)), 
the probability for a neutrino to be in ji'e) is given by 

\We\i^.,E)\'=l+^u.-^P, (46) 

where is the flavour-basis, z-direction unit vector in flavour space. Similarly, the 
probability for an antineutrino to be li?^) is given by 

mi',.E)\' = l-s^-4'^- (47) 

NFIS can also be defined from the matrix of density for a neutrino v or 
antineutrino D: 

inji Tr{g^^E^) for neutrino, 

i^n^^ Tr{(Tygi},E'^l(f) — —-^n^^ Tr{g*^ ^(t) for antineutrino, 

where 



= J J2 [Tr(^'.(P).B^O]'• (49) 

y i—x,y,z 

Here NFIS represents the average flavour of the neutrino (antineutrino) mode uj, 
and is the "net number density" of the neutrino (antineutrino) mode u; that is 
in the flavour state represented by s^. For example, if all neutrinos with energy E 
are either in pure \ve) state or in pure with number densities n^^^E and nu^,,E, 
respectively, then 

^ sgn{n^^^E - n^^,,E)^ and = l^i^^.E - n^^, ,e\- 
Using we can define the distribution function of NFIS to be 

n°J^{r) 

We emphasise that, for flavour transformation of supernova neutrinos under the single- 
angle treatment, all quantities that are proportional to neutrino number densities such 
as n°*^(r) and n^{r) are computed with the geometric factor C(r) included. Because 
we consider only forward scattering between neutrinos and background particles, /^^ 
do not change with r. We note that fi^j is usually less than the neutrino (antineutrino) 
number density with energy \5vn? I2uj\. They are equal only if all neutrinos or 
antineutrinos with this energy arc in the same flavour state. As a result, the NFIS 
distribution function /^^ is not normalised to unity but, rather, /wdw < 1. 
Under the NFIS notation equation (1) becomes 

^s^{r) = s^{r) X H^{r) = s^{r) x [wi?vac + J?matt(r-) - ^i{r) {s{r))] , (51) 
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where the misalignment between the "vacuum field" 

H^^^ = -4^ ) sin 26lv + e^P cos 20^ (52) 

-if) 

and ei indicates the mismatch between vacuum mass eigenstates and pure flavour 
states of neutrinos, the "matter field" 

i?matt(r-) = ~V2GFNAp{r)Y,{r)e^P (53) 

represents the change of refractive indices of neutrinos caused by ordinary matter, and 

{s{r)) = / dujf^s^ir) (54) 



is the average NFIS. Equation (51) shows that NFIS's are "antiferromagnetically" 
coupled with each other with strength /i(r) — 2\/2GFn1^{r). 

In vacuum, NFIS sl^ precesses around -ffvac with angular frequency |a;|. When 
projected to the -axis, this precession motion is interpreted as the oscillation of 
neutrino fiavours [see equations (46) and (47)]. 

Because of the popular usage of the neutrino polarisation vector, it shall be 
helpful to discuss the relationship between this notation and the NFIS notation. From 
equations (43) and (48) it is clear that NFIS's and neutrino flavour polarisation vectors 
are connected to each other by simple relations 

P.Mr) = 2nf{r)f^sUr) and Pp,£(r) = -2nf (r)/^sL(r). (55) 

A slightly different definition 

P^(r) cx sgn(cjW)/^sL(r) (56) 

is used in some recent literature (e.g. Raffelt and Smirnov, 2007b). With this new 
definition, Puj{r) and s^{r) are related by some constant factor (i.e. independent of 
r) which is positive for a neutrino {toSm'^ > 0) and negative for an antineutrino 
{ujSrri^ < 0). The e.o.m. for Puj{r) is 

^PUr) = M + A(r)L + ii'{r)D{r)] x P^(r), (57) 
dr 

where B = — -ffvac, ^{'r)L = — -ffmatt(''), m'(^) is different from fi{r) by a constant 
factor depending on the normalisation of Puj{r), and 

/oo 
sgn{ujSm^)P^{r) du cx (s(r)). (58) 
-oo 

We note that the orders of the cross products on the right-hand sides of equations 
(51) and (57) are different with the former closely imitating the e.o.m. for magnetic 
spins in magnetic fields. 

The notations of the polarisation vector and the NFIS are fully equivalent. 
However, there is a caveat when the "corotating frame" technique is applied (Duan 
et al, 2006c). In the absence of ordinary matter (i.e. i?matt — 0) the e.o.m. for 
the NFIS system is essentially unchanged when observed in a reference frame that 
rotates about i?vac with a constant angular frequency luq. The only difference is 
that the angular precession frequency cu of each NFIS slj(r) is shifted by —ojq. 
This is a powerful technique to analyse collective flavour transformation in neutrino 
systems. For example, consider a neutrino gas that initially consists of pure and Vr' 
with energies |5m^/3a'o| E^nd |5m^/wo|, respectively. Applying the corotating-frame 
transformation, one can see that the flavour evolution of this system is similar to that 
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of the system which initially consists of pure, mono-energetic Ve and Ve with energy 
\5m^/ujo\- We note that a corotating-frame transformation can change the sign of 
the angular precession frequency uj of polarisation vector Pi^(r), and the direction 
of Puiir) must be simultaneously inverted when this occurs. In the NFIS notation, 
however, the superficial distinction between particles and antiparticles (i.e. the positive 
and negative frequency modes) is eliminated, and s^{r) is invariant under such a 
corotating-frame transformation. In some very recent literature (e.g. Raffelt, 2008) 
polarisation vectors for antineutrinos are defined in a way similar to NFIS's and have 
their directions inverted. The caveat discussed above, therefore, disappears under this 
new definition, and ^L(f) and s^{r) are only different by a positive, constant factor 
that is proportional to /i^. 



4-.3. Adiabatic MSW-like flavour evolution 

Without knowing the results from numerical computations, one may think that the 
effects of background neutrinos can be treated simply as another potential added to 
Ve(r) in the MSW mechanism (e.g. Qian and Fuller, 1995b). In particular, if both 
p(r) and n°^(r) vary slowly with r, one would expect that a neutrino or antineutrino 
stays in the same matter eigenstate which is essentially a pure fiavour state at the 
neutrino sphere. Of course, a matter eigenstate in this case is an eigenstate of the full 
Hamiltonian which includes neutrino self-interaction. 

Similar to the energy eigenstates of an electron in the presence of a magnetic 
field, in the NFIS notation the matter eigenstates of neutrinos are represented by 
spins that are aligned or antialigned with the total effective field i?i^(r) in fiavour 
space. If supernova neutrinos experience the adiabatic MSW-like fiavour evolution 
described above, one should have 

or 

= w^AH^lttAr) - /i(r)(,sW(r))], (60) 

^L:lw--^^MO(4'Ho>, (61) 

= ^£^[c. + i/ilt,.M-M(r)(4^'(0>], (62) 

where the alignment factor e^, is constant and equals to +1 (—1) if s^{Ru) is aligned 
(antialigned) with H^{R^)^ H^{r) — \H^^{r)\, and xf^"* [i = x,y,z) stand for the 
components of vector X in the vacuum mass basis: 

4/) ^ 40 cos 26'v + e^P sin 29^ , (63) 
^;^ = ^y\ (64) 
4^) = H^^^ = -4f) sin 26'v + e^P cos 26lv . (65) 
Averaging equations (60)-(62) over all neutrino modes one obtains 

(4^H0> = i^[H^ZLA^) - m(^)(4^H0>] j J^^y (66) 
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Figure 14. Comparison of the results of the single-angle calculations shown in 
figure 12 with those in the adiabatic MSW-like and adiabatic precession solutions. 
Figure adapted from Duan et al (2007c). 



(4^H0> - / dc.-^[c. + i?il,.,(r) -^(r)(.W(r))]. (68) 
2 J-oo ti^{r) 

Equations (66) and (67) imply that, if ifmatt(?') 7^ 0, then {s'y^ (r)) — and, therefore, 
= for all neutrinos [equation (61)]. 
If neutrinos follow the adiabatic MSW-like flavour evolution, one can solve 
equations (66) and (68) for (si^^(r)) and (si^^(r)) and then solve equations (60) and 
(62) for each individual NFIS. Duan et al (2007c) have done exactly that with the 
same physical settings as those of the single-angle calculations shown in figure 12. The 

components of the average NFIS (si^^) and (s^^) = \J (s^J^)'^ + {s^^)'^ — \{s^x^)\ as 
functions of r in this solution are plotted in figure 14. For comparison the results of 
the corresponding single-angle calculations are also shown in the same figure. 

For the inverted mass hierarchy case [figures 14(a,b)] one observes that at r < 63 
km (s||^^(r)) and {s^z'\r)) in the single-angle calculation agree very well with the 
"adiabatic MSW-like solution" which is solved from equations (66) and (68). At 
r > 63 km, however, {s^^\r)) abruptly jumps out the track of the adiabatic MSW-like 
flavour evolution. For the normal mass hierarchy case [figures 14(c,d)] the neutrino 
system follows the adiabatic MSW-like solution to a larger radius. Duan et al (2007c) 
were not able to solve equations (66) and (68) beyond r ~ 91 km. Figure 14(d) 
shows that, before the neutrino system deviates from the adiabatic MSW-like flavour 
evolution, it seems to experience an MSW resonance as {s^z\r)) changes its sign at 
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r ~ 90 km. [According to equations (46) and (47) the flavour transformation of a 
neutrino is represented by the changing in the orientation of the correspond NFIS.] 

Figure 14 suggests that configuration of the neutrino system becomes unstable 
just before it departs from the adiabatic MSW-like flavour evolution. This instability 
is similar to that of a pendulum near its highest position which we will look into next. 

4- .4- Bipolar systems and flavour pendulum 

In simulating neutrino flavour transformation in the early Universe Kostclccky and 
Samuel (1993) observed an intriguing phenomenon that with an inverted mass 
hierarchy certain neutrino gases can experience "substantial flavour oscillation even 
for extremely small mixing angles" . This phenomenon can be partly understood by 
using the concept of the NFIS energy of a homogeneous, isotropic neutrino gas (Duan 
et al, 2006c): 

£{t) ^ - r diuUtosUt) ■ ^vac + ^\m)\^ (69) 

where we have assumed -ffmatt = 0- Note that the second term on the right-hand side 
of equation (69) has a positive sign because the NFIS's are "antiferromagnetically" 
coupled to each other. Using equation (51) one can show that, if and, therefore, 
\x do not vary with time t, £ must also be constant. 

As a simple example, we consider a homogeneous and isotropic neutrino gas which 
initially consists of equal numbers of pure and v^. with the same energy Eq. We 
assume that -ffmatt — and n'^J^ is constant. In the NFIS notation this system is 
represented by two NFIS's s^^o(*) and s_^^„(t) with tJo = Sm^ /2Eo. At t = NFIS's 
sl;(,(0) and s_a;o(0) are aligned and antialigned with , respectively. This is an 
example of the bipolar system which consists of two groups of NFIS's approximately 
opposing each other in direction. According to equation (69) the NFIS energy of such 
a neutrino gas, up to a constant, is 

£ = -^[^^-o(i) - ■ + ^Gpnf [sL„(i) • s^^„{t)] = const. (70) 

When n^^ is large, £ is dominated by the coupling energy between the NFIS's, 
and, because of energy conservation, s^„(t) and s^i^^^t) must remain in a bipolar 
conflguration. 

For a normal mass hierarchy ((5m^ > 0) the initial conflguration of our simple 
bipolar system becomes absolutely stable in the limit ^ (i.e. e"P ei^'). This is 
because in this limit the bipolar system is initially in the lowest energy conflguration 
where both the coupling energy between s±c^g and -ffvac and that between Sc^^ and 

5- ojg are at their minimum values. This is not the case for an inverted mass hierarchy 
((5m^ < 0) with which these two coupling energies are at their maximum and minimum 
values, respectively. Spontaneous collective neutrino oscillation is, therefore, forbidden 
to occur in the former case but is allowed in the latter. In fact, if Srri^ < and 
Gpn'^ ^ \Sm? /2Eo\, the bipolar conflguration of the NFIS's is energetically allowed 
to nearly completely flip its direction, and, therefore, both neutrinos and antineutrinos 
can almost entirely change their flavours. 

The flavour dynamics of this bipolar system is probably best illustrated using 
the pendulum analogy introduced by Hanncstad et al (2006). For this bipolar system 
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equation (51) can be recast in the form 

^ J(i) = q{t) X Mg, (71) 

J{t) = Mq{t) X ^<r(t) + am- (72) 

Equations (71) and (72) describe the motion of a pendulum with total angular 
momentum J(t) — sLo(^) + (^S-^„(t). Here, for more generality, we allow the number 
densities of antineutrinos and neutrinos to be different with a = j-u>ol fuio- The mass 
of the pendulum M = (l+a)/^ is located at position q{t) = Q{t)/Q and it experiences 
a constant gravity field g — (Qu!o/M)Hvac, where vector 

I + a ^ 

Q{t) ^ s^„{t) - as_^^{t) ^ t^o^^vac (73) 

has a constant length Q. For a ^ \ the pendulum moves like a gyroscope because it 
has a constant, non-vanishing internal spin a = J{t) ■ q(t). 

When neutrino number densities are large, the NFIS's in this example system 
maintain a bipolar configuration, and qit) ~ 2sljo(t) ~ — 2s_j^,-|(t). For a normal 
mass hierarchy {ujq > 0), g is in the same direction as that of iJvac- In the limit 
0, q(t = 0) ~ H^a,c — ff/l.9l- In other words, the flavour pendulum is near its 
lowest position. For an inverted mass hierarchy (cjq < 0), however, g is in a direction 
opposite to that of ifvac- In the limit 6^ 0, one has q{t = 0) ~ ~g/\9\, and the 
flavour pendulum is near its highest position. 

The inverted mass hierarchy case with 9^ ^ 1 is interesting. In this case, if the 
number densities of neutrinos and antineutrinos are equal, the flavour pendulum will 
always swing from near the highest position through the lowest position and back to 
its initial height. During this process bipolar system can experience significant flavour 
transformation. If the number densities of neutrinos and antineutrinos are different, 
the pendulum will undergo gyroscopic motion because of its non-vanishing spin, and 
it will not pass through the lowest position. In particular, the pendulum will remain 
at its highest position like a "sleeping top" if (e.g. Scarborough, 1958) 

> 1- (74) 

For the simple bipolar neutrino system shown above, this is equivalent to the condition 
that (Duan et al, 2007a; Hannestad et al, 2006) 

2(1 — VQ^) (-fF^o 

This condition can viewed as the division between the so-called "synchronised regime" 
and the "bipolar regime" for the following reasons. 

According to equation (51) each individual NFIS s^{t) tends to precess about 
Hvnc with its angular frequency w. At the same time, the coupling among NFIS's 
tend to make them move collectively. If the latter tendency dominates, all the NFIS's 
are locked into one block spin with angular precession frequency 

dc.VL,%^. (76) 



oo 
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This is the synchronised neutrino flavour transformation explained by Pastor et al 
(2002). A necessary condition of synchronisation is that each NFIS precesses about 
the block spin much faster than does the block spin about i/vac, i-e. 

^^m>\^sync\■ (77) 

Duan et al (2006c) pointed that a bipolar system can experience synchronised or 
bipolar (pendulum-like) oscillations depending on whether condition (77) is satisfied. 
For the simple bipolar system discussed above, the estimates given by equations (75) 
and (77) for at the boundary between the synchronised and bipolar regimes are 
different only by a constant factor. 

4.. 5. Transition to bipolar oscillations 

Pastor and Raffelt (2002) proposed that supernova neutrinos could experience 
synchronised flavour transformation when neutrino fluxes are very large. Using 
equation (77) Duan et al (2006c) determined the synchronised and bipolar regimes 
in the supernova environment. The notion of the synchronised regime has since then 
been widely adopted. However, figure 14 shows that neutrinos actually experience 
the MSW-like flavour evolution instead of the synchronised flavour transformation in 
the so-called synchronised regime when i?matt 7^ 0. In addition, in the normal mass 
hierarchy case supernova neutrinos follow the MSW-like flavour evolution well beyond 
the synchronised regime as determined in the inverted mass hierarchy case (figure 14). 

Duan et al (2007c) proposed that fiavour transformation of supernova neutrinos in 
the collective regime, i.e. neutrino self-interaction is not negligible, can be explained 
as the combination of two adiabatic solutions, i.e. the adiabatic MSW-like solution 
discussed in section 4.3 and the adiabatic precession solution to be described in section 
4.6. The fiavour pendulum model discussed in section 4.4 offers great insights into 
the question why neutrinos deviate from the adiabatic MSW-like fiavour evolution in 
supernovae. We shall elaborate on this idea in more detail. 

We note that at the neutrino sphere the most abundant neutrino species are 
and Pe, which is similar to the simplistic bipolar system represented by the fiavour 
pendulum model. In addition, as pointed out by Duan et al (2006c), the effects of 
ordinary matter can be "ignored" for collective neutrino oscillations even if the density 
of ordinary matter is large (but still low enough to be transparent to neutrinos) . The 
idea of ignoring ordinary matter is at the core of the collective flavour transformation 
of supernova neutrinos and we shall discuss it flrst. 

According to equation (51), all NFIS's tend precess to about the matter fleld 
-ffmatt with the same frequency and, therefore, -ffmatt does not break the coUectiveness 
of neutrino oscillations as 77vac does. In fact, the matter field "disappears" in the 
reference frame that rotates about iJmatt with angular frequency |i?matt|- This is 
similar to the situation where the gravity field "vanishes" in a freely-falling reference 
frame. In this corotating frame, however, iJvac is not stationary but rotates about 
—Hraatt with angular frequency l^mattl- If the matter density is large, H^^,^, the 
component of i?vac that is perpendicular to iJmatt, rotates very fast about -ffmatt- As a 
result, its effects on NFIS's average to zero if the NFIS's participate in some collective 
motion with a time scale much longer than 27r/|_ffmatt|- In other words, when the 
matter density is large, one can ignore the effects of ordinary matter for collective 
neutrino oscillations by adopting a small effective mixing angle and an effective mass- 
squared difference Sm^g — 6m^ cos20v Hannestad et al (2006) analysed the effects 
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of ordinary matter on the flavour pendulum and confirmed that matter density had 
little effect on the motion of the pendulum except for a logarithmic dependence of the 
period of this motion on the matter density. 

Let us now consider the inverted mass hierarchy case with 9^ <^ 1- We 
can estimate when bipolar oscillations may start using the analogy of the flavour 
pendulum. Near the neutrino sphere the effective total neutrino number density 
n1^{r) is large, and the mass of the pendulum M cx {n^)~^ is small. As a result, 
the pendulum with nonzero spin cr ~ (1 — a)/2 is able to "defy" gravity and stays 
at its highest position. In other words, no significant bipolar oscillation can occur 
and supernova neutrinos follow the adiabatic MSW-like flavour evolution in which 
both neutrinos and antineutrinos essentially remain in their original flavours. As 
neutrino number densities decrease with r the mass of the flavour pendulum increases. 
Eventually, the flavour pendulum becomes so heavy that its highest position is no 
longer stable [see equation (74)]. When this occurs, supernova neutrinos will break 
away from the MSW-like flavour evolution and experience the bipolar oscillation. As 
mentioned earlier, unlike in the MSW-like flavour evolution, the presence of ordinary 
matter will not affect the bipolar oscillation once it starts. 

Bipolar neutrino oscillations are represented by the gyroscopic motion of the 
flavour pendulum which is a combination of nutation and precession. The nutation 
of the flavour pendulum signifles simultaneous oscillations of s^^z, the z components 
of the NFIS's, or the oscillations neutrino survival probabilities (see flgure 12). Here 
we do not distinguish between the vacuum mass basis and the flavour basis because 
the effective mixing angle is small. The precession of the flavour pendulum stands for 
the simultaneous precession of all NFIS's about the z axis in flavour space. From the 
above discussion it is clear that in the inverted mass hierarchy case bipolar oscillations 
are insensitive to the exact value of 9^. This explains why Duan et al (2007b) obtained 
similar results for different inverted mass hierarchy schemes even when 0^ differ by 
several orders of magnitude. It is also obvious that bipolar oscillations are not sensitive 
to large matter densities in inverted mass hierarchy cases as conflrmed by Fogli et al 
(2007). 

For the normal mass hierarchy case, the flavour pendulum is initially near its 
lowest position where it is stable even when n'^{r) decreases. As a result, in the single- 
angle calculations presented in figure 12, supernova neutrinos follow the adiabatic 
MSW-like evolution until they experience significant fiavour transformation at r ~ 90 
km because of low matter density. This flavour transformation is similar to the 
resonance in the conventional MSW mechanism and is represented by flipping of the 
orientation of the NFIS's or raising of the flavour pendulum. As in the inverted 
mass hierarchy case, the flavour pendulum will undergo gyroscopic motion after 
the pendulum is raised and supernova neutrinos will, therefore, experience bipolar 
oscillations. But unlike the inverted hierarchy case, the initiation of bipolar oscillations 
depends upon the existence/efficacy of MSW-like resonances. So for the normal 
hierarchy collective neutrino flavour transformation is inevitably sensitive to the value 
of and matter densities, an expectation confirmed by Duan et al (2007b) and Fogli 
et al (2007), respectively. 

4.. 6. Adiabatic precession solution 

Comparing figures 12(c,d) and 14(b) one notices that, although individual NFIS's can 
oscillate with significant amplitudes in the z direction, [sz{r)) roughly stays constant 
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at r < 100 km. This is because, as pointed out by Hannestad et al (2006), neutrino 
lepton number 

/>oo 

L = 2 da;/^sL(r)-i/vac = 2(sW(r)) (78) 



is constant if the matter field can be "ignored" . This conservation law arises because, 
when the matter field is absent, there is only one external field Hvac in the NFIS 
system, and the e.o.m. for the NFIS's are invariant with a simultaneous rotation of 
all NFIS's about -ffvac- However, this symmetry about i?vac is usually broken because 
most configurations of NFIS systems are not invariant with rotation about Hvac- As 
a result the NFIS's tend to rotate about H^ac collectively in order to restore the 
symmetry dynamically. 

The collective mode of neutrino oscillations discussed above corresponds to the 
precession motion of a flavour pendulum. If a pendulum undergoes pure precession 
without wobbling, its motion is symmetric about the z axis although its configuration 
at any instant is not. One can imagine that such a pendulum will slowly fall but 
remain in pure precession if its weight and spin are changed infinitely slowly. Duan 
et al (2007a) studied the fiavour evolution of the simple bipolar system discussed in 
section 4.4 assuming that the fiavour pendulum stays in pure precession as n°^(r) 
decreases. They found qualitative agreement between this simple analysis and the 
numerical simulations of fiavour evolution of supernova neutrinos. 

Neutrino gases in the pure collective precession mode satisfy two conditions, 
namely, the pure precession ansatz and the adiabatic ansatz (Duan et al, 2008d). 
The precession ansatz is that at any instant all NFIS's precess about i^vac with the 
same angular frequency a;pr(n°^) which varies with n^^ , or 

^s^{r) = s^{r) X cc;p,.(nf )i?vac. (79) 
dr 

From equation (79) and the e.o.m. for NFIS's [equation (51)] one concludes that each 
NFIS slj(r) must be either aligned or antialigned with 

H^r) = H^ir) - LOp,{nf)H,,, = [cu - Wp,(nf )]i?vac - /x(nf )(s(r)). (80) 
In other words, condition 

Ur)^"f4^ (81) 

holds, where = +1 or —1 if slj(r) is aligned or antialigned with H^{r). Note that 
we have ignored i/matt in discussion of collective neutrino oscillations. 

If n1^{r) varies slowly with r, the NFIS system transits from state to state all 
of which satisfy the precession ansatz. The adiabatic ansatz is that e^^ stays constant 
during these transitions: 

= 0. (82) 

dr 

Duan et al (2006b) showed that the spectral-swap/split phenomenon (figure 13) 
can arise as a result of collective precession of NFIS's using the following simple 
argument. If neutrinos are in the pure collective precession mode, each NFIS s^^i^r) 
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must stay aligned or antialigned with H^(r) depending on the value of e^^. This is 
true even when 0. In this limit one simply have 

sl>|„rff^o = YSgn(cj - ujp^.)e^^\ (83) 

where = a;pi.(n^^ = 0). Equation (83) implies that the NFIS orientation S(j|„cff^o 
as a function of lo flips direction at uj^j.. For 0^ <C 1 this means that the neutrino 
survival probability as a function of neutrino energy jumps from to 1 or vice versa 
at energy \El\, where 

^i-^, (S4) 

The swapping point is located in the neutrino (antineutrino) sector ii El > {Ef, < 0) . 
For Ef, > the whole energy spectra of antineutrinos stay the same or get swapped 
depending on whether the neutrino mass hierarchy is normal or inverted. Raffelt and 
Smirnov (2007b) pointed out that the value of ui^^. or Ef, can be determined using 
equation (83) and the conservation of lepton number L defined in equation (78): 



L= dw/^e^sgn(w - cj|,j.). (85) 

^ — oo 

Raffelt and Smirnov (2007b) showed that, for homogeneous, isotropic neutrino 
gases or supernova neutrinos under the single-angle approximation, the pure collective 
precession mode of neutrino oscillations with any given can be solved from a small 

set of nonlinear equations. One notes that H^{r) is the total effective field for s^{r) 
in the frame that rotates about H^ac with angular frequency LO-pT:{n'^J^). Because a 
NFIS is only stationary when it is either aligned or antialigned with its total effective 
field, the precession ansatz is equivalent to the assumption that it is possible to find a 
corotating frame in which all NFIS's are stationary. In this corotating frame equation 
(81) is equivalent to 



(86) 
(87) 



where Sai.x is the projection of s^^ on e^- Here we define Cx to be the unit vector in 
the same direction of (sj_), the vector component of (s) which is perpendicular to ei^^ 
Multiplying equations (86) and (87) with and integrating them over lu one obtains 

1 =-i r diu- 



1 / , ^^Ji^OJ 



_ - dc^ , (89) 

Note that in equations (88) and (89) the minus signs arise because of the 
antiferromagnetic coupling between NFIS's. Given (L, n°*^, t^^) one can solve equations 
(78), (88) and (89) for [ujpr, (sj), {sz))- One can then solve equations (86) and (87) 
for each individual NFIS. 

Duan et al (2007c) solved the adiabatic, procession solution from equations (78) 
and (86)-(89) that corresponds to the single-angle calculations presented earlier. 
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The results are plotted in figure 14 together with the results from single-angle 
calculations and the corresponding adiabatic, MSW-like solutions. One observes 
that (sj_(r)) = |(si(r))| obtained in single-angle calculations oscillate around the 
adiabatic precession solution after neutrinos stop following the adiabatic MSW-like 
flavour evolution. The values of (sz(r)) are approximately constant when the collection 
precession mode first sets in but then changes as the densities of both ordinary matter 
and neutrinos continue to decrease with r and conventional MSW flavour conversion 
starts to intervene. 

4-7. Collective neutrino oscillations in realistic supernova environments 

We have so far summarised our understanding of collective neutrino oscillations under 
the assumptions that both the two-flavour neutrino mixing scheme and the single- 
angle treatment of a spherical supernova model are valid. The results obtained using 
these assumptions shed light on, and form the basis of, our understanding of three- 
flavour, collective neutrino oscillations in realistic supernova environments. Here we 
briefly report the progresses made in these directions. 

Duan et al (2008b) carried out the first single-angle calculations which take 
into account neutrino self-interaction and use the full three-fiavour neutrino mixing. 
These calculations are targeted at a very early epoch of an 0-Ne-Mg core-collapse 
supernova when the fast deleptonization of the PNS leads to an intense Ve burst. These 
calculations show that the neutrino energy spectra emerged from the supernova surface 
can also contain step- like features as in two-fiavour calculations. The results of these 
calculations also suggest that, like in the conventional MSW flavour transformation, 
collective neutrino oscillations in the full three-flavour mixing scheme can be factorised 
into two two-flavour mixing scenarios that occur in sequence. 

To generalise the spin-precession analogy to three flavour oscillations, Dasgupta 
and Dighe (2008) defined eight-dimensional Bloch vector , the three-fiavour version 
of the polarisation vector, as 

77°^ 1 ^ 

gc. = ^ + • A, (90) 

where {i = 1,2,..., 8) are the Gell-Mann matrices. One can define the cross 
product of eight-dimensional vectors by replacing the structure constant e^-fc of the 
SU(2) group with the structure constant fijk of SU(3) (Kim et al, 1988). The e.o.m. 
for can then be written in a form similar to equation (57). With Bloch vectors, 
three-fiavour collective neutrino oscillations can be visualised by the "ea-es" diagram, 
where (z = 1,2, ... ,8) are the unit basis vectors for the eight-dimensional fiavour 
space. Using the concept of Bloch vectors and a simplified supernova model Dasgupta 
et al (2008c) showed that the step-like features in the energy spectra of supernova 
neutrinos in Duan et al (2008b) can indeed be understood as two spectral swaps/splits 
each of which corresponds to the conservation of a neutrino lepton number. 

Independently, Duan et al (2008d) attacked the problem of three-fiavour, 
collective neutrino oscillations using a different approach. Using the correspondence 
between the operations for NFIS's and flavour matrices: 

X s^, [g^, g^,] and ■ s^' Tr(£i^g^/), (91) 

one can translate the analyses for two-flavour neutrino mixing from the spin 
representation to matrix representation. In this translation a corotating frame is 
equivalent to a unitary transformation. Similarly, the statement that a NFIS Suj is 
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stationary in such a corotating frame means that commutes with its Hamiltonian 
after an appropriate unitary transformation. Expressing the precession ansatz 
[equation (79)] and the adiabatic ansatz [equation (82)] in terms of flavour matrices and 
generahsing them to the three-flavour mixing scheme, Duan et al (2008d) showed that 
the e.o.m. for in the three-flavour scheme are invariant under independent rotations 
of the whole system about A3 and Ag when the matter fleld is absent. Two collective 
precession modes of neutrino oscillations can arise because of these symmetries which 
would naturally lead to spectral swaps/splits as in three-flavour mixing scenarios. 

While all of the above analyses on collective neutrino oscillations are based 
on the homogeneity and isotropy of the neutrino gas or the validity of the single- 
angle approximation of supernova neutrinos, a real supernova environment is clearly 
inhomogeneous and anisotropic. Some of the conclusions that are correct for 
homogeneous, isotropic neutrino gases may not be so for inhomogeneous, anisotropic 
ones. For example, Esteban-Pretel et al (2008a) recently showed that very dense 
ordinary matter can play a bigger role in suppressing collective neutrino oscillations 
in anisotropic environments than in isotropic environments. 

Raffelt and Sigl (2007) studied some simple anisotropic neutrino gases initially 
consisting of neutrinos and antineutrinos with equal numbers. They found that 
collective oscillations in such systems quickly break down and flavour equipartition 
is always achieved. We note that such a system correspond to a flavour pendulums 
with no internal spin which does not undergo precession at all. This study, therefore, 
suggests that bipolar neutrino oscillations that correspond to the nutation motion of 
the flavour pendulum are not collective in anisotropic environments. However, the 
collective neutrino oscillation mode that corresponds to the precession motion of the 
flavour pendulum arises because of a symmetry in the equations that govern flavour 
evolution of neutrinos. This symmetry also exists in inhomogeneous, anisotropic 
neutrino gases (Duan et al, 2008c). Protected by this symmetry, the collective 
precession mode of neutrino flavour transformation may not break down even in 
inhomogeneous, anisotropic environments. Indeed, Esteban-Pretel et al (2007) 
showed that single-angle-like collective neutrino oscillations appear as the neutrino- 
antineutrino symmetry is dropped which permits the precession motion of the flavour 
pendulum. 

Using the symmetry of the equations that govern neutrino flavour evolution 
Duan et al (2008c) prescribed the adiabatic, precessional flavour evolution in 
inhomogeneous, anisotropic neutrino gases, which is similar to the simpler cases with 
homogeneous and isotropic neutrino gases. The key idea is that, assuming all NFIS's 
to precess collectively about H^ac at any space-time point (instead of any instant in 
homogeneous, isotropic cases) one can flnd a corotating frame in which all the NFIS's 
at this point are stationary. Such a collective neutrino oscillation mode is described 
by a wave-like solution. As neutrinos travel from one "wavefront" to another, the 
corresponding NFIS's rotate about Hvac with a common angle which is equal to the 
phase difference between the two "wavefronts" of the solution. 

5. Summary and outlook 

While the neutrino telescopes around the globe are waiting for the neutrino shower 
from the next Galactic supernova, our rapidly-growing knowledge of neutrino flavour 
transformation in supernovae suggests that these cosmic neutrinos will reveal much 
valuable information about the supernova as well as about themselves. Through these 
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neutrinos the dynamical features of the supernova such as shocks and turbulence may 
be observed and analysed. Undetermined fundamental properties of the neutrino such 
as its mass hierarchy may be disclosed through the swapped/split neutrino energy 
spectra due to neutrino self-interaction. However, much still need to be learned about 
this subject. 

For example, collective neutrino oscillations in highly inhomogeneous, anisotropic 
environments are so far not well understood. Most of our understanding of collective 
neutrino oscillations is based on the single-angle approximation which ignores the 
anisotropic nature of the supernova environment. Meanwhile, direct numerical 
simulation of neutrino oscillations with neutrino self-interaction in such environments 
(using 2D and/or 3D matter profiles) remains to be done. The subject of turbulence 
in supernova and its impact is also uncertain. One of the imminent tasks is to combine 
the effect of the dynamic matter profile and that of neutrino self-interaction, which 
has been studied largely in parallel so far. The matter contribution is important near 
the region where G-pp{r) ~ Sm'^/{2E) and collective neutrino oscillations can occur 
if G-pnuir) > Srn^ /{2E). If the two regions are well separated, one can conduct two 
calculations in series with the results from the computation of collective oscillations as 
inputs to that on the MSW effect. This is the approximation made by Chakraborty 
et al (2008); Cava et al (2009); Kneller et al (2008); Lunardini et al (2008). In the 
case where the two regions overlap with each other, in principle one should carry out a 
comprehensive calculation with both effects included which could be a time-consuming 
task. 

In summary, we have reviewed the field of neutrino oscillations in supernovae 
focusing upon the recent discoveries of the dynamism of the pure-MSW effect and the 
emergence of neutrino self-interactions near the neutrino sphere. We have attempted 
to give those readers interested in pursuing research in this area sufficient detail that 
he or she can now begin to digest the ever-growing literature on the subject and, 
at the same time, give the casual observer a overarching understanding of these two 
phenomena. But perhaps the most general summary we can make is that the evolution 
of this field over the past few years has been both fruitful and frenetic and we fully 
expect that this trend will continue for some time yet. 
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